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We assess the performance of the Heyd-Scuseria-Ernzerhof (HSE) screened hybrid density func- 
tional scheme applied to the perovskite family LaM0 3 (M=Sc-Cu) and discuss the role of the 
mixing parameter a (which determines the fraction of exact Hartree-Fock exchange included in 
the density functional theory (DFT) exchange-correlation functional) on the structural, electronic, 
and magnetic properties. The physical complexity of this class of compounds, manifested by the 
largely varying electronic characters (band/Mott-Hubbard/charge-transfer insulators and metals), 
magnetic orderings, structural distortions (cooperative Jahn- Teller like instabilities), as well as by 
the strong competition between localization/delocalization effects associated with the gradual filling 
of the t2 g and e g orbitals, symbolize a critical and challenging case for theory. Our results indicates 
that HSE is able to provide a consistent picture of the complex physical scenario encountered across 
the LaM0 3 series and significantly improve the standard DFT description. The only exceptions 
are the correlated paramagnetic metals LaNi03 and LaCuOs, which are found to be treated better 
within DFT. By fitting the ground state properties with respect to a we have constructed a set of 
'optimum' values of a from LaSc0 3 to LaCu0 3 : it is found that the 'optimum' mixing parameter 
decreases with increasing filling of the d manifold (LaSc0 3 : 0.25; LaTiOs & LaV0 3 : 0.10-0.15; 
LaCr0 3 , LaMn0 3 , and LaFe0 3 : 0.15; LaCo0 3 : 0.05; LaNi0 3 & LaCu0 3 : 0). This trend can be 
nicely correlated with the modulation of the screening and dielectric properties across the LaM0 3 
series, thus providing a physical justification to the empirical fitting procedure. 

PACS numbers: 71.27.+a, 71.15.-m, 71.10.-w, 71.28.+d, 71.20.Bc, 78.30.Er, 75.25.-j 



I. INTRODUCTION 

The physics of transition metal perovskites with gen- 
eral chemical formula ABO3 (where A is a large cation, 
similar in size to O 2- and B is a small transition metal 
(TM) cation) has attracted and challenged the interest 
and curiosity of the material science community for many 
decades, due to huge variety of complex phenomena aris- 
ing from the subtle coupling between structural, elec- 
tronic and magnetic degrees of freedom. The high de- 
gree of chemical flexibility and the localized (i.e. not 
spatially homogeneous) character of the dominant TM 
partially filled d states lead to the coexistence of several 
physical interactions - spin, charge, lattice and orbital 
- which are all simultaneously active. The occurrence 
of strong lattice-electron, electron-spin and spin-orbit 
couplings causes several fascinating phenomena, includ- 
ing metal-insulator transitions 1 * 2 -, superconductivity^, 
colossal magnetoresistance 4 ^, multiferroicity^, bandgaps 
spanning the visible and ultraviolet^, and surface chem- 
ical reactivity from active to inert^. When the addi- 
tional degrees of freedom afforded by the combinatorial 
assemblage of perovskite building blocks in superlattices, 
heterointerfaces and thin films are introduced the range 
of properties increases all the more, as demonstrated by 
the recent several remarkable discoveries in the field of 
oxide heterostructures 1 ^. Tunability and control of these 
intermingled effects can be further achieved by means of 
external stimuli such as dopin g 11 : 12 , pressur e) 13 : 14 , tem- 
perature and magnetic or electric field s 15 : 16 , thereby en- 



hancing the tailoring capability of perovskites for a wide 
range of functionalities. This rich array of behaviors 
uniquely suit perovskites for novel solutions in different 
sectors of modern technology (optoelectronics, spintron- 
ics, piezoelectric devices and (photo) catalysis), for which 
conventional semiconductors cannot be usecU^r— . 

Theoretical studies of TM perovskites, aiming to de- 
scribe and understand the underlying physical mecha- 
nisms determining their complex electronic structures 
have been mainly developed within two historically dis- 
tinct solid state communities, model Hamiltonian s 21 : 22 
and first principles^, which in the recent years have ini- 
tiated to fruitfully cross-connect each other methodolo- 
gies towards more general schemes such as DFT+DMFT 
(Density Functional Theory^ 4 . + Dynamical Mean Field 
Theory^ 5 - - — ), with the aim to overtake the individual 
limitations and to improve the applicability and pre- 
dictive power of electronic structure theor y 28 ' 29 . Model 
Hamiltonians approaches adopt simplified lattice fermion 
models, typically the celebrated Hubbard model, inspired 
by the seminal works of Anderson^, Hubbard 3 ^ and 
Kanamori 3 -^ in which the many-body problem is solved 
using a small number of relevant bands and short-ranged 
electron interactions. These effective models can solve 
the many-body problem very accurately, also includ- 
ing ordering and quantum fluctuations, but critically 
depends on an large number of adjustable parameters 
(which can be in principle derivable by first principles 
schemes 3 - 3 . - — ), and its applicability is restricted to finite- 
size systems 2 -^ 9 -. In DFT the intractable many-body 
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problem of interacting electrons is mapped into a sim- 
plified problem of non-interacting electrons moving in an 
effective potential throughout the Kohn-Sahm scheme^, 
and electron exchange-correlation (XC) effects are ac- 
counted by the XC potential which is approximated us- 
ing XC functional such as the Local Density Approxi- 
mation (LDA), the Generalized Gradient Approximation 
(GGA) et similari(S§-. As the name suggests, in DFT 
the ground state properties are obtained only from the 
charge density, and this makes DFT fundamentally differ- 
ent from wavefunction-based approaches as the Hartee- 
Fock method, the simplest approximation to the many- 
body problem which includes the exact exchange but no 
correlational. Though DFT has been widely and success- 
fully used in the last 40 years in solid-state physics and 
quantum-chemistry to calculate structural data, energet- 
ics and, to a lesser extent, electronic and magnetic prop- 
erties, it suffers of fundamental difficulties mostly due 
to the approximate treatment of XC effects. This draw- 
back is particularly severe when DFT is applied to the 
so called strongly correlated systems (SCSs), whose pro- 
totypical examples are transition metal oxides (TMOs). 
A systematic improvement of these XC-rclated deficien- 
cies in DFT is essentially impossible, but several 'beyond- 
DFT' approaches have been proposed which deliver much 
more satisfying results. The most renewed ones are the 
DFT+U 38 ; SIG22ri2, hybrid functional^, and GWil 

In this article we applied the screened hybrid func- 
tional introduced by Hcyd, Scuseria, and Ernzerhof 
(HSE)i£ to study the structural, electronic and magnetic 
properties of the series of 3d TMO perovskites La7\/03, 
with M ranging from Sc to Cu. This is a rather challeng- 
ing family of compounds for electronic structure methods 
for several reason o 35 ' 46 " — (see Table |TJ (i) it encompasses 
band, Mott-Hubbard (MH) and charge-transfer (CT) in- 
sulators as well as correlated metals (the last two mem- 
bers of the series: LaNiOs and LaCuOa)^; (ii) differ- 
ent type of antiferromagnetic (AFM) ordcrings are en- 
countered across the series (A-type, C-type and G-type, 
graphically represented in Figfl]), but also non-magnetic 
(NM, LaSc0 3 ) and paramagnetic (PM, La(Co->Cu)0 3 
systems 4 ^; (iii) the dominating electronic character varies 
from d° to d 8 , and ranges from ti g je g localization (with 
variable crystal filed splitting between ti g and e g states) 
to more spatially dclocalized d-orbitals 4 ^; (iv) the crys- 
tal symmetry spans orthorhombic (O), monoclinic (M), 
rhombohcdral (R) and tetragonal (sketches of the crystal 
structures is given in Figf2]) characterized by a different 
level of structural distortions (Jahn- Teller (JT: staggered 
disproportionation of the M-0 bondlengths), GdFe03- 
type (GFO: collective tiltings and rotations of the oxygen 
octahedra), monoclinic angle (3). 

Before describing the method and presenting the result 
we briefly recall previous ab initio investigations of this 
set of compounds performed using conventional DFT and 
"beyond-DFT" methodologies. The most widely stud- 
ied member of this family is certainly the classical JT- 
GFO distorted Mutt-Hubbard AFM insulator LaMn0 3 , 
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FIG. 1. (Color online) Schematic representation of the typical 
magnetic orderings for the perovskites. 
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FIG. 2. (Color online) The structures of perovskite oxides 
studied in this paper. P„ mn for LaSc03, LaTiOa, LaCrO;?, 



LaMn03, and LaFeOa, P 
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for LaV03, i?3 C for LaCoOs 



and LaNiOs, and P^/m for LaCuOs, respectively. The large 
(green), medium-sized (blue), and small (red) balls denote La, 
M, and O atoms, respectively. Plot done using the VESTA 
visualization program^. 



but also other compounds have received significant at- 
tention, in particular LaTiOs and LaV03, and to a lesser 
extent, LaFcOs, LaCo03, LaNiOs, and LaCu03. Rela- 
tively scarce studies on the band insulator LaSc0 3 are 
present in literature. 

^32^46^49^53-56,61,68,76,79,82 . The sem i na l works of 

the Terakura group in the late 90 a 48 ' 49 ' 53 " — have exten- 
sively assess the performance of LDA for the LaA/03 
series (M=Ti-Cu), and revealed that LDA is unable to 
predict the observed insulating ground state for the first 
members (LaTiOs and LaVOs), wrongly favors a non 
magnetic solution for LaTi03 and severely underesti- 
mate the insulating gap in LaCrOs, LaMn03, LaFcOs, 
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TABLE I. Summary of the fundamental ground state properties of ha,MOs: (i) Crystal structure: 0=orthorhombic, 
M=monoclinic, R=rhombohedral, and T=tetragonal; (ii) Transition metal (TM) spin-projected electronic configuration and 
(line below) corresponding oxidation state, (iii) electronic character: I=insulator and M=metal; Magnetic ordering: NM=non- 
magnetic, different type of AFM arrangements (see FigfT} , and P=paramagnetic. 
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and LaCo03. The situation does not improve using the 
GGA££. However, the recent GGA-based re-exploration 
of the electronic properties of LaCr03 by Ong et alp- 
has reported that a good agreement with experiment can 
be achieved, upon a proper (re) interpretation of the op- 
tical spectra. It should be noticed that all these results 
were obtained using the experimental geometries. The 
very few structural optimizations at DFT level, mostly 
focused on LaMn03, have shown that though LDA/GGA 
reproduce the experimental volume within 1-3 % 64 i 76 i 82 , 
the lattice distortions associated with the JT and GFO 
instabilities are significantly underestimated. For com- 
pounds with more dclocalized 3d electrons like LaNiOs 
the LDA performance get better as recently reported by 
Guo et al2^. 

DFT+ [ ^35. 46.48. 51. 52. 54. 55. 58. 59. 62-64. 66-68. 70-74. 76-79. 84 . 

In some cases the drawbacks of LDA and GGA in 
treating localized partially filled d states can be adjusted 
by introducing a strong Hartree-Fock like intra-atomic 
interaction U properly balanced by the so-called double- 
counting (dc) correction. The resulting LDA(GGA)+U 
energy functional can be written as2£ 

E to t(n, n) = -Edft (n) + Eu F (h) - E dc (n) (1) 

where n is the operator for the number of electrons occu- 
pying a particular site and n is its expectation value. This 
expression can be written in terms of the direct (U) and 
exchange (J) contributions, which lead to a set of slightly 
different LDA(GGA)+U energy functional depending 
on the way the dc term is constructed^ 5 .. Among the 
numerous applications of DFT+U to LaAf O3 the study 
of I. Solovyev and coworkers represents the most com- 
prehensive and systematic one^. There it is found that 
LDA+U conveys a substantially improved description of 
the bandstructure of LaM03 from LaTi03 to LaCu03 
with respect to conventional DFT, though the results 
critically depend on the specific treatment of localiza- 
tion effects in the 3d manifold. By applying U-correction 
to t- 2g electrons only the authors show that LaTiO-3 and 
LaVC>3 are correctly predicted to be insulating, thus cur- 
ing the deficient LDA picture. At variance with DFT 
LaTiOs is found to be magnetic, but with a magnetic 
moment twice larger than the experimental one. The 



bandgap of early (LaTi03, LaVOs) and late (LaCoOa) 
LaM03 members which have a predominant ti g charac- 
ter are better described than the e g compounds LaMnO,-? 
and LaFeOs, for which an on site U applied to the entire 
3d is needed to improve the agreement with experiment. 
The values of the gap clearly depends on the value of 
the U parameter, as discussed by Z. Yang et al£Z. By 
fitting the U using the measured gap as reference quan- 
tity, these authors have shown that the best agreement 
with experiment is achieved for U progressively increas- 
ing from 5 cV (LaCr0 3 ) to 7 cV (LaNi0 3 ), about 2 eV 
smaller than those computed by I. Solovyev using con- 
strained LDA^ 8 -. Similarly to the standard LDA case, few 
attempts have been made to optimize the structural pa- 
rameters at DFT+U \eve&£L2&IL22.-. (i) LaTi0 3 . H.S. 
Ahm^ 7 - and coworkers have shown that the application of 
LDA+U (U=3.2 eV) systematically increase the (under- 
estimated) LDA lattice parameters of LaTi03 and the in- 
ternal distortions, thus improving the overall agreement 
with experiment, (ii) LaMn0 3 . using the Perdew-Burkc- 
Ernzerhof^ (PBE) approximation with an on-site effec- 
tive U=2 cV T. Hashimoto et a/.— have performed a full 
(volume and internal coordinates) structural optimiza- 
tion in LaMnOs and demonstrated that, unlike GGA, 
GGA+U accounts well for the experimental JT and tilt- 
ing distortions; (iii) LaCo0 3 . H. Hsu et alJQ and A. 
Laref et al^Z- have shown that LDA+U describes well the 
lattice parameter, rhombohedral angle and atomic coor- 
dinates of LaCoOs; better agreement with experiment is 
obtained using a self-consistent U2^ rather than a fixed 
U value of «7-8 eVZI. (iv) LaNi0 3 . The work of G. 
Guo et alZ^ on LaNiOs reported that for this correlated 
metal LDA+U (U=6 eV) delivers geometrical data very 
similar to the already satisfactory LDA ones (though, as 
already pointed out, LDA does a better job in predicting 
the electronic properties). 

ff j 160 ! 65 ' 87 : The application of a purely Hartee-Fock 
(HF) procedure, i.e. including an exact treatment of the 
exchange interaction and neglecting electron correlation, 
has been extensively investigated by T. Mizokawa and A. 
Fujimori^ and by I. Solovyev^ 5 -. Though the HF method 
suffers for the absence of electron correlation which is re- 
flected by its tendency to overestimate the magnitude of 
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band gaps (which can be cured by including the correla- 
tion effects beyond the HF approximation), these stud- 
ies show that HF can qualitatively explain the ground 
state electronic and magnetic properties of this class of 
magnetic oxides. Important exceptions are LaNiOs and 
LaCuC>3, which are found to be FM insulator (LaNiOs) 
and G-type AFM insulator (LaCuOs), in contrast with 
the observed PM metallic ground state. Another crit- 
ical case for HF and in general for electronic structure 
methods is the origin of the type-G AFM ordering in 
LaTiOj^^^: in Ref|H9 the authors report that the 
stabilization of the G-type arrangement can be achieved 
by fixing the Ti — O — Ti angle to approximately the 
experimental value. The resulting magnetic moment, 
downsized by spin orbit interaction effect, results in good 
agreement with the measured value, but the calculated 
bandgap is dramatically wrong, about 2.7 cV, against 
the measured value of 0.1 eV— . The results of Refroa go 
to the opposite direction: the magnetic ground state re- 
mains wrong even upon inclusion of correlation effects, 
but the bandgap, 0.6 eV, is in much better agreement 
with experiment. A similar trend is also observed for 
LaV0 3 . 

Hybrid Functional£Z£L£&l&12.-£L22.: An alternative 
methodology to DFT and HF which has attracted a 
considerable attention in the solid state physics and 
chemistry communities in the last two decades is the 
so called hybrid functional approach. Originally intro- 
duced by A.D.J. Becke in 1993^ 3 -, the hybrid functional 
scheme relays on a suitable mixing between HF and 
local / semilocal (LDA/GGA) DFT theories, in which a 
portion of the exact non — local HF exchange 
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is mixed with the complementary LDA/GGA 
local / semilocal approximated exchange E^ FT (r). 
The resulting general hybrid XC kernel E^ y ^ rld (decom- 
posed over its exchange (X) and correlation (C) terms) 
can be written in the form: 
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where the mixing factor a controls the amount of exact 
E^ F incorporated in the hybrid functional. Similarly to 
DFT+U (which makes use of the HF-like intra-atomic in- 
teraction U, as recalled above) hybrid functionals tends 
to correct the LDA/GGA dclocalization error and to pro- 
vide a better description of TMO with partially filled d 
and / states. The advantages with respect to DFT+U 
is that hybrid functionals (i) do not suffer from the dou- 
ble counting term (see Eq. Q]) and, even most impor- 
tantly, (ii) use an orbital dependent functional acting on 
all states, extended as well as localized (in the DFT+U 
method, the improved treatment of exchange effects is 
limited to states localized inside the atomic spheres, and 
usually limited to the partially filled TM shell). Though 



both schemes problematically depends on semiempirical 
parameters such as U and J in DFT+U and the mixing 
factor a in hybrid functionals, many attempts have been 
made to overcome these difficulties^—. 

Though sparse in literature, hybrid functionals 
studies of LaA/03 are increasing in the last few 
years^^^^^^^-iffi. Applications of the 
renowned Becke, 3-parameter, Lee- Yang-Parr B3LYP 
functional^ 3 , to LaMnO^^^ have shown that this 
method properly favors the type-A AFM ground state 
and provides an accurate description of the band gap, 
magnetic coupling constants, and Gibbs formation 
energies. The only structural optimization of the 
JT-distorted structure, however, deliver lattices con- 
stants which deviate by 5% from experiment. We 
have recently reported that HSE performs very well 
in predicting the ground state properties of LaMn03, 
including the optimized structural parameters, and that 
the data are slightly dependent on the actual value of the 
mixing facto r 35 i 82 . D. Gryaznov et al. have successfully 
studied the structural and phonon properties of LaCo03 
using the PBEO (Perdew-Ernzerhof-Burke)i2i hybrid 
functional and reported a substantial improvement with 
respect to conventional DFT. The application of HSE 
and PBEO functionals to LaNi03, conversely, turned 
out to give poor agreement with the experimental 
photocmission spectroscopy (PES); this is in line with 
precedent unsatisfactory HSE/PBEO results obtained 
for other itinerant magnetic metals^^. The influence of 
the non-local exchange on the electronic properties of 
LaTi0 3 has been investigated recently by F. Iori and 
coworkersiS^. By adopting the experimental structure 
these authors clarified that the improved description 
of HSE over DFT+U is due to a correct re-positioning 
of the O p states, and show that by fixing the mixing 
parameter a to its "standard" value, 0.25, the bandgap 
and the magnetic moment are significantly overestimated 
with respect to measurements. 

ff/ C 103 ' 104 : Another approach to correct the self- 
interaction (SI) LDA/GGA problem is the self- 
interaction correction method 3 -^— , in which an ap- 
proximated (atomic-like and orbitally averaged) self- 
interaction is subtracted from the LDA XC functional. 
Though conceptually different from LDA+U (in LDA+U 
an additional effective Coulomb term is added to the 
LDA/GGA functional), the SIC method is often prag- 
matically viewed as a generalized LDA+U approach in 
which the atomic SI plays the role of the U-iS 3 -. Several 
implementations of the SIC scheme have been proposed, 
characterized by a different level of complexity in treating 
the SI term and from the different underlying computa- 
tional framework 3 ^—, but all demonstrated an apprecia- 
ble accuracy in predicting and interpreting the electronic 
structure of a vast range of systems, including SCSs and 
TMOs±2£ 

A valid illustration of the performance of the SIC 
method is supplied by the results obtained for LaTi03 re- 
cently discussed by Filippctti et alJ^. By assuming the 
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experimental cell parameters SIC finds the correct AFM 
type-G insulating ordering and delivers internal struc- 
tural distortions close to the experimental ones. As a 
downside, however, the magnitudes of the bandgap (1.6 
cV) and magnetic moment (0.89 /ie) are substantially 
larger than the corresponding measured values («0.2 eV 
and ~0.5 /ib, respectively). Other SIC applications to 
the LaMC>3 series are limited, to our knowledge, to the 
ideal undistorted cubic phase of LaMnOy^, for which a 
stringent comparison with experiment is difficult to do. 

GH ^ 35 ' 69 ' 71 : We finally recall the main achievements on 
LaAf O3 acquired using the GW approximation, a compu- 
tational method fundamentally different from both DFT 
and HF. GW is configured to reflect and to treat the 
quasi-particle nature of electrons on the basis of Green's 
function many-body perturbation theory^, by explicitly 
accounting for the non-local and frequency-dependent 
self-energy (E) in a suitably rewritten Schrddingcr-likc 
equation. In the GW approximation E is approximated 
to the lowest order term of the Hedin's equation, and can 
be written as: 



E = iGW (4) 



where G is the Green's function and W is the dynam- 
ically screened Coulomb kernel. In the most widely 
used single-shot Go Wo approximation both G and W 
arc treated in an unperturbed manner, but with in- 
creasing computer power self-consistent or partially self- 
consistent GW schemes are becoming more and more 
possibl o 35 ' 106 ' 107 . Due to the extensive computing time 
required to perform GW-like calculations, only few GW 
data are available in literature for complex systems. 
Among these, the works of Y. Nohara et aZi 69 i 71 repre- 
sent a very comprehensive example of a systematic ap- 
plication of GW to LaM03 starting from preconvergent 
LDA+U wavef unctions. These authors have obtained ex- 
cellent agreement with experimental spectra, but proba- 
bly due to the uncertainties connected to the choice of U 
in preparing the initial wavefuncitons the values of the 
computed bandgaps deviate significantly from the ex- 
perimental estimations, especially for LaTi03, LaV03, 
and LaCo03. Good agreement with experiment has 
been also obtained for LaMn03 using a partially self- 
consistent GWo approach, in this case starting from 
GGA wavefunction^ 5 .. 

The paper is organized as it follows. In Section |TT] 
we illustrate the computational method and its techni- 
cal aspects; in Section IIIII we report the results on the 
structural optimization (Section IIII A[) and electronic & 
magnetic properties (Section IIII B|) . A more general dis- 
cussion on the observed trends and behaviors is devel- 
oped in Section HVl and finally in Section [V] we draw our 
summary and conclusions. 



II. COMPUTATIONAL ASPECTS 

All calculations were performed using the Vienna ab 
initio Simulation Package (VASP) 108 i 109 employing DFT 
and hybrid-DFT approaches within the projector aug- 
mented wave metho d 110 ' 111 and the PBE paramctriza- 
tion scheme^ for the XC functional. In the screened 
hybrid-DFT HSE approach adopted in this study, part of 
the short-range (sr) PBE exchange functional is replaced 
by an equal portion of exact HF exchange, according to 
the general prescription: 

E« S C E = aE BF ^+(l-a)E P / E - sr ^+E P / E ^+E EBE 

(5) 

where (i, controls the range separation between the sr and 
long-range (lr) part of the Coulomb kernel (1/r, with 
r = \r — r'|), decomposed over long (L) and short (S) 
terms: 

- = S^r) + Lp(r) = — - — + (6) 

The reason to include a screening parameter (i is moti- 
vated by the computational effort required in computing 
the spatial decay of the HF X interaction. In the refined 
HSE06 hybrid functional, fj, is set equal to 0.20 A -1 which 
corresponds to the distance 2//i at which the HF X inter- 
actions starts to become negligible. For /x=0 the PBEO 
functional is recovered^, whereas for /i — s- 00 HSE be- 
comes identical to PBE. Beside the computational cost, 
the main beneficial consequence of the inclusion of a 
screening strategy in PBEO is that screened hybrids can 
give access to the metallic state, which is unaffordable 
by unscreened PBEO-likc hybrids, and generally improves 
the accuracy of bandgaps for semiconductors^^. 

The mixing parameter a, determining the amount of 
exact non-local HF X included in the hybrid XC func- 
tional is usually set to 0.25^. In this case the PBE func- 
tional is recovered for a=0. 

Thus, the HSE06 depends by construction on two pa- 
rameters, \i and a. Though their standard values are rou- 
tinely used in solid state calculations, it is to be expected 
that they may vary from material to materia l 97 ' 113 , or 
that they may be property-dependen t 57 ' 1 14 i 1 15 . Unfor- 
tunately, a rigorous first principle procedure to deter- 
mine the choice of these parameters does not exist. The 
conventional value a~l/A is determined by perturba- 
tion theoryiSi. The choice /j, = 0.20 A^ 1 has proven to 
be a practical compromise between computational cost 
and quality of the results^. Considering that most of 
the tests and fitting procedures have been performed 
taking as a reference atomic or molecular energetical 
and structural properties, the direct acquisition of these 
standard values in extended solid state system is not 
straightforward 101 ' 116 . 

By linking hybrid density functional theory with many- 
electron XC self-energy E within a GW framework it 
has been proposed that the mixing factor a can be 
interpreted as the inverse of the dielectric constant 
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Co d 113 ' 117 Based on this idea, an approximated recipe 
to determine the optimum value of a can be obtained: 



a op t ~ — (7) 



which depends solely on the dielectric constant and on 
the "unknown" factor of proportionality. It follows 
straightforwardly that for metal (e^ = oo) a opt is equal 
to zero. Several limitations affect this practical rule and 
degrade its ab initio natur o 119 , above all an accurate cal- 
culation of the dielectric constant, which is presently very 
difficult in particular for complex TMOs. 

Following this line of thought other strategies have 
been introduced to overcome this problem invoking 
density-functional estimators'^ in the spirit of the Tran 
and Blaha functionali^i, which furnishes parametric ex- 
pressions inevitably dependent on the specific mate- 
rial dataset, usually limited to monoatomic and binary 
semiconductors^ 7 .. To complicate the situation even fur- 
ther, there is some amount of arbitrariness in transfer- 
ring the a op t ~ — relation from unscreened PBEO-likc 
hybrids to screened ones like HSE, where screening is al- 
ready present in some form in the range separation con- 
trolled by the screening factor \x. These complications 
become particularly cumbersome when one moves from 
'standard' monoatomic and binary semiconductors to the 
more complex ternary TM perovskites. As a matter of 
fact, due to the absence of a systematic study on the in- 
fluence of a in this class of compounds, the large majority 
of hybrid functionals studies on ternary TMOs have been 
performed using the standard 1/4 compromise, though 
there are neither fundamental nor practical justifications 
for this choice. 

Thus, in order to shed some light on the role of a in 
a representative class of ternary TMOs with a largely 
varying degree of screening and competition between lo- 
calization/delocalization effects we have performed our 
HSE calculations using 4 different values of a: (i) low 
mixing (strong screening): 0.10 (HSE-10), 0.15 (HSE- 
15), (ii) standard mixing: 0.25 (HSE-25), and (iii) high 
mixing (low screening) : 0.35 (HSE-35). The careful anal- 
ysis of structural, electronic and magnetic properties will 
allow us to draw some general trends which should serve 
as a guidance for future HSE applications. 

Technical setup: The plane wave cutoff energy was 
set to 300 eV. A 4x4x4, 6x6x6, 8x8x8 Monkhorst- 
Pack fc-points grids were used to sample the Brillouin 
zones for ¥ nrna /¥ 2l /b, Bg c , an d P4/m structures, respec- 
tively. Structural optimization was achieved by relaxing 
the volume, lattice parameters, lattice angles and inter- 
nal atomic positions throughout the minimization of the 
stress tensor and forces using standard convergence cri- 
teria. 



III. RESULTS 

This section is subdivided into two parts which are de- 
voted to the presentation of the structural (Sec. IIII Al) . 
and electronic and magnetic (Sec. IIII B[) properties, re- 
spectively. In each subsection we will summarize the spe- 
cific results obtained for each member of the LaM03 se- 
ries, and in the next section (Sec. IIV[) we will provide a 
more reasoned discussion on the general trends observed 
across the series. 

As anticipated in the introduction hybrid function- 
als can be simplistically viewed as an orbital-dependent 
DFT+U approach in which the on-site electron-electron 
interaction parameter U is replaced by the parametric in- 
clusion of a portion of the exact HF exchange quantified 
by the mixing factor a. In DFT+U calculations, the U 
is usually either tuned to fit some specific physical prop- 
erty (i.e. bandgap, magnetic moment, volume, etc.), or 
calculated within constrained-LDA procedures 9 -^— . In 
contrast, most of the available HSE- based calculations 
present in literature are done at fixed mixing parame- 
ter a=0.25. This might erroneously convey the idea of a 
minor role played by the mixing factor or, even more fun- 
damentally misleading, that HSE is a purely ab initio (i.e. 
parameter free) scheme. As already discussed previously, 
in the last few years the modeling community has started 
to address this issu o 57 ' 97 ' 113 but the amount of available 
data are still very limited, in particular for complex oxide. 
It is therefore instructive to briefly recall a few results on 
the choice of the U in DFT+U studies of transition metal 
oxides in order to possibly formulate some expectations 
on the behavior of the mixing parameter a in HSE. A 
good example to start with are transition metal monox- 
ides (TMOs: MnO, FeO, CoO and NiO), where the TM 
possesses the oxidation state 2+ (M 2+ ). Several LDA+U 
studies have shown that a U between 6 and 8 eV can pro- 
vide an accurate enough prediction of band gaps for all 
TMO o 38 ! 122 . Going from M 2+ to M 3+ the number of the 
localized electrons decreases. Thus, it might be expected 
that the magnitude of the Coulomb interaction increases 
due to the contraction of the spatial extension of the of 
the 3c? (M 3+ ) wave functions'^ 3 - However, by compar- 
ing M 2+ and LaM 3+ 03 photoemission data it can be 
unambiguously concluded that the effective Coulomb in- 
teraction decreases in M 3+ compounds^^ 3 --^. 

Under the assumption that in LaM 3+ 03 the t^g elec- 
trons are localized and the e g electrons are itinerant, 
Solovyevi£ has explained this apparent contradiction by 
invoking the strong screening associated with the e g elec- 
trons. Indeed, the computed value of U for the t 2g shell 
in LaM 3+ 03 are significantly reduced with respect to U 
for the d states in M 2+ 0. The strength of the screening 
depends on the filling of the e g orbital: it is strong at half- 
filling and less efficient when the e g are nearly empty of 
occupied. The results of Solovyev indicate that this ^g-U 
approach reproduces sufficiently well the main features of 
early (Ti-V-Cr) and late (Co-Ni) LaM03 compounds but 
fails for LaFeOs (much too small bandgap and magnetic 
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moment) and LaMnOs (small gap). Clearly, effects other 
than e g itinerancy contribute to the strength of the U, 
such as the screening from non-3d electrons, M (3d / 4s)- 
0(2p) hybridization, and lattice relaxation which can ex- 
plain the discrepancy between self-consistent +U meth- 
ods and experiments. 

A fitting-U approach can selectively adjust the com- 
parison with the experimental gap (not for LaCr0 3 ) at 
the expense of a rigorous description of the position of 
the e g , tig and 0(p) sub-bands (i.e. the 'correct' value of 
the bandgap can arise from an fundamentally incorrect 
artificial electronic structure). This failure prevents any 
physically sound specification/understanding of the (MH 
or CT) character of the gap: in RefJH, for instance, the 
gap of LaMnC>3 is found to be predominantly CT like, 
in discrepancy with the actual situation (LaMnOs is a 
MH insulator with a gap opened between occupied and 
empty e g states, partially hybridized with O p states). 
Furthermore, the 'optimum' U's resulting from fitting- 
U schemes do not seem to reflect the observed M 2+ to 
M 3+ U-reduction, which is an additional sign of the in- 
adequacy of such a procedure. 

Considering that standard HSE (a=0.25) performs 
quite well for TM monoxide a 112 ' 129 we can expect that 
a smaller value will turn out to be more appropriate 
for reproducing the ground state properties of La-A/03. 
Furthermore, given the full-orbital character of HSE we 
may expect that hybridization effects and screening from 
non 3d electrons will be better described as compared 
to DFT+U. Finally, we should point out that the choice 
to perform a complete structural optimization at each 
considered value of a allows for a more genuine account 
of the structural contribution to the screening which is 
disguised in frozen lattice (atomic positions fixed to ex- 
perimental ones) calculations. 



A. Structural Properties 

As already mentioned in the introduction four dif- 
ferent crystal symmetries are encountered across the 
LaM03 series (see Fig. (i) Orthorhombic P nma 

for LaSc0 3 , LaTi0 3 , LaCr0 3 , LaMn0 3 , and LaFe0 3 ; 
(ii) Monoclinic P2 1 /b for LaVOs; (iii) Rhombohcdral i?g c 
for LaCoOs and LaNiOs; and (iv) Tetragonal P 4 / m for 
LaCu03. All these different structures share the same 
octahedral perovskitic building block MOq, character- 
ized by one central TM metal surrounded by two apical 
(Oi) oxygen atoms and four planar (O2) oxygen atoms. 
Depending on the specific compound, the MOq octa- 
hedra can undergo two kinds of structural distortions: 
the JT distortion, manifested by a short (s) and long (1) 
M-O2 in-plane distances and medium (m) M-Oi verti- 
cal ones (along the octahedral axis), and the GFO tilt- 
ing of the M-OT- M and M - O^- M 180° angles. 
The cooperative JT distortion is usually measured in 
terms of the JT modes Q2=2(Z-s)/-\/2 and Q3=2(2m-Z- 
s)/-\/6. In our full structural relaxation we have Opti- 
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FIG. 3. Experimental trend of Volume (V), average tilting 
angle (8), and JT distortions (Q2 and Q3), for the 
series from M=Sc to M=Cu. The corresponding trend of the 
tolerance factor (t=(RA+Ro)/v / 2(RM+Ro), and Rm is also 
shown. 



mized the volume (V), lattice parameters a, b, and c, the 
monoclinic/rhombohcdral angle j3, as well as all inter- 
nal atomic coordinates (this clearly includes all relevant 
GFO and JT structural parameters M-Ch-M (6»x), 
M-O^-M (9 2 ), Q 2 , and Q 3 ). 

A graphical summary of the observed trend of the 
most relevant structural parameters is given in Fig. [3] 
The progressive reduction of the Volume from Sc to Cu 
is clearly associated with the almost monotonically de- 
crease of the M ionic radius Rm, whose size is deter- 
mined by the competition between the size of the 4s shell 
(where extra protons arc pulled in) and the additional 
screening due to the increasing number of 3c? electrons: 
adding protons should lead to a decreased atom size, but 
this effect is hindered by repulsion of the 3d and, to a 
lesser extent, 4s electrons. The V/Rm curves show a 
plateau at about half filling (Cr-Mn-Fe) indicating that 
for this trio of elements these two effects are essentially 
balanced and atom size does not change much. The vol- 
ume contraction's associatedjvith a rectification of the 
average (M - Oi - M+M - 2 - M)/2 tilting angle 9, 
which follows very well the evolution of the tolerance fac- 
tor <=(Ra+Ro)/v / 2(Rm+Ro) (where Ra, Rm and Ro 
indicate the ionic radius for La, M=Sc-Cu and O, respec- 
tively) . This indicates that the tolerance factor is indeed 
a good measure of the overall stability and degree of dis- 
tortion of pcrovskitc compounds. Clearly, the value of t 
is well within the range of stability set to 0.78 < t < 1.05. 
The bottom panel of Fig. [3] conveys the message that Q2 
and Q3 assume non negligible values for LaMnOs only, 
confirming that JT distortions are predominant in per- 
ovskitcs containing cations such as Cu 2+ and Mn 3+ in 
their octahedral cation site. 

In the following subsections we will report on the full 
structural optimization of LaM03 at PBE and HSE (for 
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TABLE II. Structural data for LaScOg. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing factor) and the avail- 
able (room temperature) experimental data taken from Ref. 
Il3dl . The relative error (in brackets, in %) and the mean 
absolute relative error (MARE, %) is also supplied. 



Expt. HSE-35 HSE-25 HSE-15 HSE-10 PBE 



v (A 3 ) 


266.09 


262.02 


263.48 


265.12 


265.99 


267.90 






(1.5) 


(1.0) 


(0.4) 


(0.0) 


(0.7) 


a (A) 


5.787 


5.764 


5.780 


5.794 


5.798 


5.810 






(0.4) 


(0.1) 


(0.1) 


(0.2) 


(0.4) 


6(A) 


8.098 


8.050 


8.061 


8.076 


8.088 


8.108 






(0.6) 


(0.5) 


(0.3) 


(0.1) 


(0.1) 


C (A) 


5.678 


5.647 


5.655 


5.666 


5.672 


5.686 






(0.5) 


(0.4) 


(0.2) 


(0.1) 


(0.1) 


Sc-o m (A) 


2.104 


2.091 


2.096 


2.100 


2.103 


2.108 






(0.6) 


(0.4) 


(0.2) 


(0.0) 


(0.2) 


Sc-o, (A) 


2.140 


2.095 


2.101 


2.108 


2.109 


2.115 






(2.1) 


(1.8) 


(1.5) 


(1.4) 


(1.2) 


Sc-o s (A) 


2.096 


2.082 


2.086 


2.091 


2.093 


2.098 






(0.7) 


(0.5) 


(0.2) 


(0.1) 


(0.1) 


fi n 


148.39 


148.42 


148.08 


148.14 


148.19 


148.18 






(0.0) 


(0.2) 


(0.2) 


(0.1) 


(0.1) 


02 n 


146.29 


149.98 


149.95 


149.55 


149.68 


149.53 






(2.5) 


(2.5) 


(2.2) 


(2.3) 


(2.2) 


MARE 




1.0 


0.8 


0.6 


0.5 


0.6 


Q 2 


0.063 


0.018 


0.021 


0.024 


0.023 


0.023 


Q 3 


-0.023 


0.004 


0.005 


0.000 


0.003 


0.002 



different values of a) and will provide a one-to-one com- 
parison with available experimental data, also in terms 
of the mean absolute relative error (MARE, not given for 
the very small quantities Q2 and Q3). 

1. LaScOs 

LaScOa crystallizes with a P nmQ orthorhombic struc- 
ture, and shows the largest tilting instabilities of all 
LaM0 3 series (147.3°)±2a. The JT parameter Q 2 and 
Q3 are almost zero (0.063 and -0.023, respectively) and 
as a consequence the Sc-0 bondlcngth disproportionation 
is negligible: both planar and vertical Sc-0 bondlcngths 
are all rj 2.1 A. The computed structural data are col- 
lected in Table |nj All methods deliver a quite satisfactory 
description with an overall MARE less than 1%. PBE 
supplies the best agreement with measurements (MARE 
= 0.5 %). The most critical quantities for theory are Sc- 
O; and 62, for which relative errors larger than 1% and 
2% are found, respectively. We can thus conclude that 
for an accurate account of the structural properties of 
LaScOa it is not required to apply bcyond-DFT meth- 
ods. As we will see in the next section, this is not the 



TABLE III. Structural data for LaTiOa. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing fa ctor) and the avail- 
able experimental data taken from Ref. 1 13 ll . The relative 
error (in brackets, in %) and the mean absolute relative error 
(MARE, %) is also supplied. 



Expt. HSE-35 HSE-25 HSE-15 HSE-10 PBE 



v (A 3 ) 


249.17 


250.03 


250.00 


250.98 


251.17 


250.61 






(0.3) 


(0.3) 


(0.7) 


(0.8) 


(0.6) 


a (A) 


5.589 


5.599 


5.597 


5.612 


5.617 


5.646 






(0.2) 


(0.1) 


(0.4) 


(0.5) 


(1.0) 


fc(A) 


7.901 


7.931 


7.909 


7.915 


7.913 


7.929 






(0.4) 


(0.1) 


(0.2) 


(0.2) 


(0.4) 


C (A) 


5.643 


5.631 


5.648 


5.651 


5.650 


5.598 






(0.2) 


(0.1) 


(o.i) 


(o.i) 


(0.8) 


Ti-O s (A) 


2.028 


2.038 


2.034 


2.033 


2.029 


2.018 






(0.5) 


(0.3) 


(0.2) 


(o.o) 


(0.5) 


Ti-O, (A) 


2.053 


2.059 


2.069 


2.071 


2.068 


2.033 






(0.3) 


(0.8) 


(0.9) 


(0.7) 


(1.0) 


Ti-O m (A) 


2.032 


2.027 


2.023 


2.028 


2.030 


2.029 






(0.2) 


(0.4) 


(0.2) 


(o.i) 


(0.1) 


Oi (°) 


153.78 


153.12 


152.98 


153.54 


154.25 


158.47 






(0.4) 


(0.5) 


(0.2) 


(0.3) 


(3.0) 


02 n 


152.93 


152.64 


152.54 


152.59 


152.97 


156.37 






(0.2) 


(0.3) 


(0.2) 


(o.o) 


(2.2) 


MARE 




0.31 


0.33 


0.35 


0.31 


1.07 


Q 2 


0.029 


0.046 


0.065 


0.062 


0.054 


0.006 


Q 3 


-0.023 


-0.008 


-0.021 


-0.027 


-0.031 


-0.021 



case for the electronic properties. 



2. LaTiOa 

Similarly to LaScC>3, the low temperature space group 
of LaTi0 3 is P Ilrna with small JT distortions due to the 
low JT activity of the single t\ g orbital and large GFO 
distortions caused by the large size difference between Ti 
and La ions . 131 ' 132 Though the overall PBE MARE is 
only 1%, the relaxed structure parameters given by HSE 
functionals arc appreciably better (MARE w 0.3%) than 
PBE, regardless the amount of exact HF exchange, as 
summarized in Table Mil The PBE errors mostly arise 
from an incorrect description of the tilting angles, which 
are by far (« 3%) overestimated with respect to the low- 
temperature experimental data^ 3 ^. As for the volume, 
PBE furnishes a nice optimized value, which is improved 
by going to large a HSE setups (HSE-35 HSE-25). 
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TABLE IV. Structural data for LaV03. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing factor) and the avail- 
able (T=10 K) experimental data taken from Ref. 1 1331 . The 
relative error (in brackets, in %) and the mean absolute rela- 
tive error (MARE, %) is also supplied. 



Expt. HSE-35 HSE-25 HSE-15 HSE-10 PBE 



v (A 3 ) 


241.10 


240.31 


241.33 


242.20 


242.45 


241.64 






(0.3) 


(0.1) 


(0.5) 


(0.6) 


(0.2) 


a (A) 


5.5917 


5.562 


5.582 


5.622 


5.637 


5.613 






(0.5) 


(0.2) 


(0.5) 


(0.8) 


(0.4) 


b(A) 


7.7516 


7.801 


7.787 


7.729 


7.713 


7.729 






(0.6) 


(0.5) 


(0.3) 


(0.5) 


(0.3) 


o(A) 


5.5623 


5.538 


5.552 


5.574 


5.577 


5.570 






(0.4) 


(0.2) 


(0.2) 


(0.3) 


(0.1) 


n°) 


90.13 


89.93 


90.16 


90.16 


90.18 


90.03 






(0.2) 


(0.0) 


(0.0) 


(0.1) 


(0.1) 


v 2 -o s (A) 


1.979 


2.019 


1.993 


1.974 


1.966 


1.962 






(2.0) 


(0.7) 


(0.3) 


(0.7) 


(0.9) 


V2-0, (A) 


2.039 


2.007 


2.019 


2.055 


2.054 


2.012 






(1.6) 


(1.0) 


(0.8) 


(0.7) 


(1.3) 


v 2 -o s (A) 


1.979 


2.004 


1.999 


1.984 


1.991 


2.011 






(1.3) 


(1.0) 


(0.3) 


(0.6) 


(1.6) 


Vi-o s (A) 


1.978 


1.969 


1.989 


1.972 


1.965 


1.961 






(0.5) 


(0.6) 


(0.3) 


(0.7) 


(0.9) 


Vi-O; (A) 


2.042 


2.007 


2.026 


2.063 


2.059 


2.013 






(1.7) 


(0.8) 


(1.0) 


(0.8) 


(1.4) 


Vi-O m (A) 


1.989 


1.997 


1.997 


1.982 


1.990 


2.013 






en d\ 

(U.I) 


en d\ 


en d\ 

(V.I) 


(n 1 1 

(U.1J 


Cl 9 s ! 
(l.Z) 


0i n 


156.74 


155.83 


155.88 


156.65 


157.70 


160.15 






(0.6) 


(0.5) 


(0.1) 


(0.6) 


(2.2) 




157.83 


156.23 


156.90 


157.05 


157.06 


158.76 






(1.0) 


(0.6) 


(0.5) 


(0.5) 


(0.6) 


02i n 


156.12 


157.08 


156.23 


156.28 


156.51 


158.26 






(0.6) 


(0.1) 


(0.1) 


(0.2) 


(1.4) 


MARE 




0.82 


0.48 


0.35 


0.51 


0.98 


Q21 


0.085 


0.003 


0.028 


0.101 


0.088 


0.001 


Q31 


-0.050 


0.023 


-0.027 


-0.075 


-0.093 


-0.080 


Q22 


0.074 


0.015 


0.042 


0.114 


0.098 


0.002 


Q32 


-0.060 


-0.054 


-0.037 


-0.082 


-0.097 


-0.084 



3. LaVOi 

LaVC>3 is the only member of the Lail/03 scries dis- 
playing a monoclinic structure with ¥ 2l /b space groupie. 
The unit cell contains two inequivalent V sites (Vi and 
V2), which sit in the center of GFO distorted octahe- 
dra not subjected to significant JT distortions. Due to 
the occurrence of two different V atoms in the unit cell 
two different sets of V-0 bondlcngths and tilting angle 
62 (O21 and 622) arc identified. The comparison between 



TABLE V. Structural data for LaCr03. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing factor) and the avail- 
able (room temperature) experimental data ta ken fro m Re f. 
1 1 3-41 (similar structural data can be found in Refs ll35l and ll36T ). 
The relative error (in brackets, in %) and the mean absolute 
relative error (MARE, %) is also supplied. 



Expt. HSE-35 HSE-25 HSE-15 HSE-10 PBE 



V (A 3 ) 


235.02 


233.45 


234.74 


236.13 


236.69 


237.70 






(0.7) 


(0.1) 


(0.5) 


(0.7) 


/ 1 1 \ 

(1.1) 


a (A) 


5.483 


5.478 


5.509 


5.494 


5.531 


5.512 






(0.1) 


(0.5) 


(0.2) 


(0.9) 


(0.5) 


b (A) 


7.765 


7.752 


7.766 


7.776 


7.785 


7.795 






(0.2) 


(0.0) 


(0.1) 


(0-3) 


(0-4) 


c (A) 


5.520 


5.498 


5.487 


5.527 


5.531 


5.533 






(0.4) 


(0.6) 


(0.1) 


(0.2) 


(0.2) 


Cr-Oi (A) 


1.977 


1.973 


1.977 


1.982 


1.984 


1.983 






(0.2) 


(0.0) 


(0.3) 


(0.4) 


(0.3) 


Cr-O m (A) 


1.972 


1.971 


1.975 


1.979 


1.981 


1.985 






(o.i) 


(0.2) 


(0.4) 


(0.5) 


(0.7) 


Cr-O s (A) 


1.970 


1.970 


1.975 


1.979 


1.980 


1.984 






(0.0) 


(0.3) 


(0.5) 


(0.5) 


(0.7) 


0i (°) 


158.14 


158.29 


158.10 


159.76 


157.71 


158.51 






(o.i) 


(o.o) 


(1.0) 


(0.3) 


(0.2) 


02 (°) 


161.32 


159.72 


159.60 


157.66 


159.73 


159.30 






(1.0) 


(1.1) 


(2.3) 


(1.0) 


(1.3) 


MARE 




0.30 


0.30 


0.60 


0.51 


0.61 


Q2 


0.003 


0.001 


0.001 


0.000 


0.001 


0.001 


Qa 


0.010 


0.004 


0.004 


0.004 


0.005 


-0.002 



the low-temperature experimental data and the theoret- 
ical values are collected in Table IIV1 The general sit- 
uation is similar to LaTi0 3 : the PBE MARE, 1%, is 
about twice larger than the average HSE MARE. The 
PBE relative errors are large for the tilting angles and 
V-0 bondlengths, but rather small for volume and lattice 
constants. HSE leads to slightly better data, in partic- 
ular in the range 0.10 < a < 0.25, but the volume and 
lattice constants are reproduced less accurately than at 
PBE level. 



4. LaCrOz 

The structural data for P„ ma LaCr03 are collected in 
Table [V] The full three-fold degenerate t2 g shell inhibits 
completely any tendency to JT distortions but the size 
difference between La and Cr drives a substantial GFO- 
like tilting of the CrOe octahedra. Also in this case PBE 
performs as well as HSE-10 and HSE-15 (MARE = 0.6 
%). The overall MARE is further reduced to 0.3 % for 
larger values of a (HSE-25 and HSE-35). 
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TABLE VI. Structural data for LaMnO;?. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing factor) and the avail- 
able (T=4.2 K) experimental data taken from Ref. 1 1371 . The 
relative error (in brackets, in %) and the mean absolute rela- 
tive error (MARE, %) is also supplied. 





Expt. 


HSE-35 HSE-25 HSE-15 HSE-10 


PBE 


V ( k 3 \ 

V (A ) 


0/1 Q KT 

z4o.o / 


243.98 


245.82 


247.36 


244.24 


O/1/l 91 
Z44.Z1 






(0.2) 


(0.9) 


(1.6) 


(0.3) 


CO 31 


n (k\ 

a (A) 


o.ooz 


5.526 


5.537 


5.553 


5.661 


k Kan 
O.Ooy 






(0.1) 


(0.1) 


(0.4) 


(2.3) 


CO 71 


h ( k\ 
b (A) 


K 7 AO 
0. ( 4Z 


5.789 


5.817 


5.820 


5.594 


K fi97 
O.OZ / 






(0.8) 


(1.3) 


(1.4) 


(2.6) 


(2 01 


r. (k\ 

c (A) 




7.628 


7.633 


7.653 


7.712 


7 7QQ 






(0.5) 


(0.5) 


(0.2) 


(0.6) 


(1.0) 


ivin-w m (A) 


1 0^7 

i .yo t 


1.954 


1.957 


1.962 


1.979 


1 QQ9 

i.yyz 






(0.2) 


(0.0) 


(0.3) 


(1.1) 


CI 81 


Mn-Oi (A) 


2.184 


2.204 


2.214 


2.213 


2.134 


2.063 






(0.9) 


(1.4) 


(1.3) 


(2.3) 


(5.5) 


Mn-O s (A) 


1.903 


1.899 


1.905 


1.914 


1.923 


1.971 






(0.2) 


(0.1) 


(0.6) 


(1.1) 


(3.6) 


0i (°) 


154.3 


154.78 


154.35 


154.36 


153.96 


155.85 






(0.3) 


(0.0) 


(0.0) 


(0.2) 


(1.0) 




156.7 


154.38 


154.08 


154.17 


157.59 


157.71 






(1.5) 


(1.7) 


(1.6) 


(0.6) 


(0.6) 


MARE 




0.52 


0.66 


0.81 


1.22 


1.90 


Q 2 


0.398 


0.431 


0.437 


0.423 


0.298 


0.131 


Q3 


-0.142 


-0.159 


-0.167 


-0.165 


-0.080 


-0.041 






5. 


LaMnOz 









LaMnOs is the most critical case across the LaA/03 
series, due to the concomitant occurrence of both GFO 
and JT structural distortions, the latter originating by 
the intrinsic instabilities associated with the orbital de- 
generacy in the e g channel of the Mn 3+ cation. The 
lattice constants and atom positions of 0-P nma LaMnOs 
were fully optimized at both PBE and HSE level within 
an AFM-A magnetic ordering, though, as we will dis- 
cuss in the next section, PBE is not able to catch the 
correct magnetic ground state and rather favors an FM 
arrangement. The results are listed in Tabic fVTl In this 
case PBE docs not supply a satisfactory account of the 
structural properties, reflected by a quite large MARE 
(« 2%) significantly larger than the corresponding HSE- 
10 (1.22), HSE-15 (0.81), HSE-25 (0.66), and HSE-35 
(0.52). The major obstacle for PBE is the correct pre- 
diction of the JT distortions: (i) the relative error for 
the M-0 bondlength disproportionation is as high as 5.5 
%, and (ii) Q2 and Q3 are found to be one third of the 
measured values. The serious underestimation of Q2 and 
Q3 at PBE level has important consequences on the clcc- 



TABLE VII. Structural data for LaFcO.i. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing factor) and the avail- 
able (room temperature) experimental data taken from Ref. 
Il38l . The relative error (in brackets, in %) and the mean 
absolute relative error (MARE, %) is also supplied. 





Expt. 


HSK-35 




HSE-15 


HSE-10 


PBE 


V (A ) 


o/io on 


z4U.oy 


9/19 fl8 

z4z.Uo 


9/i/i ri9 

Z44.UZ 


9/1 ^ no 
z4o.uy 


9/1£ /17 
Z4D.4 ( 






C1 01 


CO 31 


CO 5~1 


CO Q) 


CI 51 


n ( k\ 

a (A) 


O.ODO 


O.OoU 


o.ooy 


O.OO / 


K ^^7 

0.00 / 


O.Olo 






(0 61 


Co 11 


CO 41 


CO 1 1 


CI 01 


h CA1 
b (A) 


i .oD4 


7 S9Q 

1 .ozy 


7 8/19 
( .o4Z 




( .0D0 


7 878 
1 .0 I 






CO 31 


CO 21 


CO P> 

(U.L) 


CO 2) 


CO 31 


r- ( k\ 

c (A) 


K ^7 


O.OOo 


0.04O 


O.OOU 


O.DUO 


O.Ouo 






(ft 1 "\ 


(ft X\ 


(ft ft\ 


(ft Q\ 


(ft 9> 


Vo ci, ( A 1 
re-\Ji (A) 


z.uuy 


9 nm 


9 ftft« 
Z.UUD 


9 019 
Z.U1Z 


9 m ^ 

Z.U10 


9 018 






(ft 41 

(V-'i) 


(ft 11 


CO 1} 


CD ^"i 


CO 41 


Fe-Oi (A) 


2.009 


2.002 


2.010 


2.017 


2.024 


2.032 






(0.3) 


(0.0) 


(0.4) 


(0.7) 


(1.1) 


Fe-O s (A) 


2.002 


1.995 


2.002 


2.008 


2.013 


2.018 






(0.3) 


(0.0) 


(0.3) 


(0.5) 


(0.8) 


fi n 


155.66 


155.95 


155.52 


155.16 


155.07 


154.83 






(0.2) 


(o.i) 


(0.3) 


(0.4) 


(0.5) 


02 n 


157.26 


157.10 


156.69 


156.29 


155.78 


155.21 






(0.1) 


(0.4) 


(0.6) 


(0.9) 


(1.3) 


MARE 




0.38 


0.16 


0.30 


0.56 


0.79 


Q 2 


0.010 


0.010 


0.011 


0.013 


0.016 


0.020 


Qs 


0.005 


0.004 


0.000 


-0.001 


-0.006 


-0.011 



tronic properties; we will discuss this issue in the next 
section. We can anticipate that the deficient treatment of 
Q2 and Q3 prevents the opening of the band gap thereby 
leading to a metallic solution. HSE-10 improves the es- 
timations of Q2 and Q3 with respect to PBE, and with 
increasing a the MARE get progressively reduced down 
to 0.52 for cv=0.35. The inaccuracy of local functional in 
reproducing the JT distortions was recently overviewed 
by Hashimoto et a/.—. In particular, these authors have 
pointed out that DFT+U can only supply a semiquan- 
titative account of JT changes if the structure is fully 
relaxed (including volume). This is also valid for purely 
unrestricted HF approaches. All other non-JT related 
quantities are equally well described by both methodolo- 
gies, with relative error generally smaller than 1 %, apart 
from the tilting angles which suffers of slightly larger de- 
viations (1-1.5 %). 

6. LaFe0 3 

The crystal of LaFeC>3 is orthorhombic with P nma 
space groupie. In the high-spin Fe 3+ configuration 
tz g *^eg^ the JT distortions are completely suppressed. 
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TABLE VIII. Structural data for LaCo03. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing factor) an d the avail- 
able (T=4 K) experimental data taken fr om Ref. Il40l (room 
temperature data can be found in Ref. 1 1411 ). The relative 
error (in brackets, in %) and the mean absolute relative er- 
ror (MARE, %) is also supplied. Here 0i=O - Co - O, and 
6>2=Co — O — Co. For LaCo03 we have optimized the struc- 
ture using a reduced 0.05 mixing parameter and obtained the 
following data: V = 111.87 (A 3 ) (1.5), a = 5.354 (A) (0.2), 
Co-Oi = 1.940 (A) (0.8), O-Co-O = 88.07 (°) (0.5), and Co- 
O-Co = 161.00 (°) (1.3); MARE = 0.82. 





Expt. 


HSE-35 


HSE-25 


HSE- 15 


HSE- 10 


PBE 


V(A 3 ) 


110.17 


107.78 


109.02 


110.39 


111.09 


114.11 






(2.2) 


(1.0) 


(0.2) 


(0.8) 


(3.6) 


a (A) 


5.342 


5.314 


5.328 


5.343 


5.348 


5.405 






(0.5) 


(0.3) 


(0.0) 


(0.1) 


(1.2) 


P (°) 


60.99 


60.70 


60.87 


61.05 


61.20 


60.99 






(0.5) 


(0.2) 


(0.1) 


(0.3) 


(0.0) 


Co-o (A) 


1.924 


1.904 


1.915 


1.925 


1.932 


1.948 






(1.0) 


(0.5) 


(0.1) 


(0.4) 


(1.2) 


fi n 


88.55 


88.96 


88.72 


88.49 


88.27 


88.50 






(0.5) 


(0.2) 


(0.1) 


(0.3) 


(0.1) 


02 n 


163.08 


165.56 


164.06 


162.88 


161.82 


162.45 






(1.5) 


(0.6) 


(0.1) 


(0.8) 


(0.4) 


MARE 




0.95 


0.42 


0.09 


0.44 


0.92 



TABLE X. Structural data for LaCuOg. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing factor) and avai lable 
(room temperature) experimental data taken from Ref. 1 1441 . 
The relative error (in brackets, in %) and the mean absolute 
relative error (MARE, %) is also supplied. 







noiL-oO 


noijj-zo 


no-Cj-io 




r DUi 


V (A ) 




ou.u / 


ou.oo 




^7 98 


^7 






(3 2) 


(2 7) 


Cl 5) 


Cl 1) 


CO 21 


a (A) 


3.819 


3.821 


3.832 


3.844 


3.850 


3.867 






(0.1) 


(0.3) 


(0.7) 


(0.8) 


(1.3) 


c(A) 


3.973 


3.840 


3.840 


3.861 


3.865 


3.869 






(3.3) 


(3.3) 


(2.8) 


(2.7) 


(2.6) 


Cu-o ; (A) 


1.986 


1.920 


1.920 


1.930 


1.933 


1.934 






(3.3) 


(3.3) 


(2.8) 


(2.7) 


(2.6) 


Cu-o s (A) 


1.909 


1.911 


1.916 


1.922 


1.925 


1.934 






(0.1) 


(0.4) 


(0.7) 


(0.8) 


(1.3) 


MARE 




2.01 


2.01 


1.70 


1.64 


1.59 


Q 2 














Q3 


0.125 


0.016 


0.007 


0.014 


0.013 


0.002 



The optimized structural data, collected in Table IVII| 
show that PBE overestimates the volume by 1.5 % but 
describes well all other parameters, leading to a relatively 
small MARE of 0.79 %. HSE predicts a better volume, 
especially for a equal to 0.15 and 0.25, but in general 
HSE improves only marginally the PBE results. 



TABLE IX. Structural data for LaNiOg. Comparison be- 
tween the optimized parameters calculated using PBE and 
HSE (with different values of the mixing factor) an d the avail- 
able (T=1.5 K) experimental data taken from Ref. Il42i . Here 
0i =0 - Ni~- O, and 2 =Ni -~CT- Ni. The relative error (in 
brackets, in %) and the mean absolute relative error (MARE, 
%) is also supplied. 





Expt. 


HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 


v (A 3 ) 


112.48 


112.02 


112.47 


113.42 


113.83 


115.20 






(0.4) 


(0.0) 


(0.8) 


(1.2) 


(2.4) 


a (A) 


5.384 


5.377 


5.380 


5.393 


5.392 


5.415 






(0.1) 


(0.1) 


(0.2) 


(0.1) 


(0.6) 


P (°) 


60.86 


60.85 


60.95 


61.01 


60.21 


61.16 






(0.0) 


(0.1) 


(0.2) 


(1.1) 


(0.5) 


Ni-0 (A) 


1.933 


1.930 


1.935 


1.941 


1.947 


1.953 






(0.2) 


(0.1) 


(0.4) 


(0.7) 


(1.0) 


oi n 


88.78 


88.78 


88.63 


88.55 


88.30 


88.41 






(0.0) 


(0.2) 


(0.3) 


(0.5) 


(0.4) 


92 n 


64.82 


164.79 


163.77 


163.43 


162.12 


163.02 






(0.0) 


(0.6) 


(0.8) 


(1.6) 


(1.1) 


MARE 




0.10 


0.19 


0.43 


0.84 


0.92 



7. LaCo0 3 

At low-temperature, LaCoOs possesses a slightly 
GFO-distorted pcrovskite structure with rhombohcdral 
symmetry (i?g c space group) 139 ' 140 , characterized by a 
rhombohcdral angle of 60.99° (see Fig. [5]). The struc- 
tural data are given in Table IVIHl The best agreement 
with experiment is achieved by HSE-15, but also HSE-10 
and HSE-25 lead to relative errors < 1%. The HSE- 
25 set of data is in good agreement with previous PBE0 
results^ PBE performs not bad (MARE below 1%), but 
overestimates too much the volume (+3.5% with respect 
to experiment). 



8. LaNiOs 

LaNiOa crystallizes with a rhombohedral structure 
with moderate GFO-like distortions^. The fully opti- 
mized structural parameters are listed in Table IPX! Sim- 
ilarly to the isostructural LaCoOa, also in this case PBE 
gives a large volume (+2.4%) but all other structural 
quantities are well reproduced (our data are in line with 
the previous calculation of Guo et al^>). Within HSE, 
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larger the amount of HF exchange is included, the more 
the MARE is reduced: from 0.84 % (HSE-10) down to 
0.1 % (HSE-35). 

9. LaCuOa 

LaCuC>3 is the only member of the LaMC>3 family dis- 
playing a tetragonal structure (P4/ m ), which can be suit- 
ably tuned to a rhombohedral one under different oxygen 
pressure conditions In this paper we only examine the 
tetragonal phase. The small elongation of the Cu06 oc- 
tahedron associated with the tetragonal form induces a 
local JT-type distortion, manifested by four equatorial 
Cu-0 bonds close to 1.909 A and two apical bonds to 
1.986 A. The relaxed structure parameters are shown in 
Table IXl From the structural data it is clear that LaCuC>3 
represents the most challenging compound of the whole 
series for both level of theory, with MARE well above 1%. 
PBE provides the overall best agreement with experiment 
(MARE = 1.59 %) but produces an almost cubic struc- 
ture, dissimilar from the observed tetragonal one. Hybrid 
functionals opens up a small structural disproportiona- 
tion between long and short Cu-0 bondlcngths which is 
however insufficient to stabilize a well defined tetragonal 
form: the lattice parameter c is still very badly accounted 
for (relative error of about 3%). 

10. Concluding remarks 

Summing up the results presented in this section we 
can draw the following conclusions: in general the struc- 
tural properties of the LaM03 series are sufficiently well 
described by standard PBE, i.e. with overall MARE 
smaller than 1 %. The only two exceptions are (i) 
LaMn.03: HSE is essential to treat correctly the JT dis- 
tortions which are a crucial ingredient to find and explain 
the A-AFM ordered insulating orbitally ordered state of 
this compound, (ii) LaCuC>3: neither PBE nor HSE 
are capable to deliver MARE smaller than 1.5%. The 
amount of non-local HF exchange does not have a deci- 
sive and univocal effect on the structural properties apart 
for the d° band insulator LaScC>3 for which the results 
get worse with increasing a. In all other cases, an im- 
provement over PBE results is obtained for all values of 
a tested in the present study, and the standard 0.25 com- 
promise seems to appear a reasonable choice. However, 
as we will discuss in the next section with this value of a 
the band gaps are found to be exceedingly overestimate, 
with respect to the measured ones. 

B. Electronic Structures and Magnetic Properties 

The focus of the present section is the presentation of 
the electronic (density of states (DOS), bandstructures 
and band gaps) and magnetic (magnetic moment m and 



Experimental Structure Optimized Structure 




Energy (eV) Energy (eV) 

FIG. 4. (Color online) ^-projected DOS of non-magnetic 
LaScOa with experimental (left) and relaxed (right) structure 
based on PBE and HSE (HSE-35, HSE-25, HSE-15, HSE-10) 
functionals. The shadow area indicates the total DOS. 



magnetic energies for different spin orderings) results for 
the entire LaM03 series, given for both the experimental 
and the fully optimized structures. Here we are particu- 
larly interested on the modifications induced in the calcu- 
lated quantities by the value of the mixing parameter a, 
from (PBE) to 0.35 (HSE-35). To this aim, following 
the outline adopted in the previous section we will se- 
quentially discuss the electronic and magnetic structure 
case by case. A more general discussion on the evolution 
of the chemical and physical properties of LaMOs per- 
ovskites from M=Sc to M=Cu will be provided in the 
next section. 



1. d°: LaSc0 3 

LaSc03 is a non magnetic band insulator with the d° 
(Sc 3+ ) electronic configuration, and an optically mea- 
sured band gap of about 6.0 eV opened between the 0-2p 
valence band and the Sc 3d unoccupied band.— Our 
calculations confirm this picture as seen from the density 
of states shown in Fig. [1J but the band gap value pre- 
dicted at PBE (3.81 eV and 3.92 eV for the experimen- 
tal and fully optimized structures, respectively, in agree- 
ment with previous calculations^) seriously underesti- 
mates the experimental value. The HSE data collected 
in Table IXII indicate that the correct value of the gap is 
recovered by admixing 25% of HF exchange. Clearly, the 
band gap increases with increasing HF percentage, but 
the DOS (see Fig. 0| always provide the same qualitative 
O-p/Sc-d picture. The band dispersion associated with 
the 25% choice given in Fig. [5] shows that the band gap 
is direct and located at T, but given the flatness of the 
topmost occupied bands (O-p) and, to a lesser extent, 
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Experimental Structure 



Optimized Structure 



TABLE XI. The band gap A (eV) of LaSc0 3 calculated 
within PBE and HSE (HSE-0.10, HSE-15, HSE-25, and HSE- 
35) using both the Experimental and the relaxed structures 
(see Table ITT)) . Other theoretical values are also listed for com- 
parison, along with the experimental measurements. 



Theory 




Optimized Structure 




HSE-35 


HSE-25 HSE-15 HSE-10 


PBE 


6.495 


5.730 4.995 4.635 


3.915 




Experimental Structure 




HSE-35 


HSE-25 HSE-15 HSE-10 


PBE 


6.435 


5.685 4.920 4.560 


3.810 




Other works 




LDA 






3.98 a 






Experiment 


~6.0\ 5.7 C 



a RefM, 'Ref.i, c Ref.ii£. 
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FIG. 6. (Color online) /-projected DOS of AFM-G ordered 
LaTiOs with experimental (left) and relaxed (right) structure 
based on PBE and HSE (HSE-35, HSE-25, HSE-15, HSE-10) 
functionals. The shadow area indicates the total DOS. 
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FIG. 7. Bandstructure of LaTiOs computed at HSE level 
(a=0.15) using the optimized structure. 



FIG. 5. Bandstructure of LaScOa computed at HSE level 
(a=0.25) using the optimized structure. 



the Sc-d bands at about 6 eV, the value of the (direct) 
band gap does not change much in the entire k-space. 
This is in agreement with the experimental optical spec- 
tra which show a sudden and very intense onset of the 
optical conductivity at 6 eV— . 



2. d 1 : LaTiOs 



LaTiC>3 is a G-AFM MH insulator with a magnetic 
moment of about 0.5 /j, e 131]14g , in which the single 3d 
electron occupies one Ti t2 g orbital. The physics of 
the orbital degree of freedom has attracted considerable 
attentio n 131 i 147 . This nominal t\ g configuration gives rise 
to a distinctive orbitally-ordered ground state character- 
ized by a very small band gap of 0.1-0.2 eV— «^ which 
has spurred a lot of theoretical study aiming to clarify the 
physics underlying this peculiar behavio r 65 i 103 i 149 ~ i££. 
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TABLE XII. The band gap A (eV), magnetic moment m 
(/is/Ti), magnetic energy (given with respect to the FM en- 
ergy, in meV) of LaTi0 3 . calculated by PBE and HSE (HSE- 
35, HSE-25, HSE-15, HSE-10) using both the experimental 
and relaxed structures (Tab. IIII|I . The gaps in brackets are 
for the G-type which is not the most favorable ordering for 
a=0.10. Other theoretical values are also listed for compari- 
son, along with the experimental measurements. 



Theory 
Optimized Structure 





HSE-35 


HSE-25 


HSE-15 


HSE-10 
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2.835 
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-17 






Experimental Structure 
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A 
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-36 


-53 


-70 


-49 


C-AFM 


-35 


-30 


-7 


32 


-21 


G-AFM 


-63 


-68 


-65 


-52 


-5 






Other works 








LDA 


LDA+U 


GW 


HF 




Band gap 




0.49° , 0.5 6 


1.00° 


2.7 e ,0.6 / 








0.4 c , 0.16 d 








m 




0.68°, 0.92 6 


0.68 a 


0.55 e ,0.76 / 





0.52 9 , 0.7' 
0.58 d 



Experiment 


A 




0.1 h c, 0.2 l m 


m 




0.45 ! a, 0.57 m 6 


ORefcS, 6 Ref4£ 


c Ref.-£2, 


d Ref.£L, e Ref.^, ^Ref.-^, 3 ReiM, 


"Ref.1, 'Refiifi, 


'Ref.iSi, 


m Ref.I2I 



In agreement with previous theoretical findings we find 
that local DFT, though it furnishes a very good descrip- 
tion of the structural properties, is incapable to repro- 
duce the MH insulating state and wrongly stabilize an 
AFM-A magnetic ordering. HSE delivers a coherent pic- 
ture which is however a dependent as summarized in 
Table IXIII and Figj6] Regardless the value of the mix- 
ing parameter, HSE predicts an insulating ground state. 
For a=0.10 HSE conveys a band gap of about 0.1/0.2 
eV (depending on whether the experimental or the fully 
optimized structure is adopted), in excellent agreement 
with experiment. However, we found that for a < 0.10 
HSE finds the AFM-A as the most favorable magnetic 



ordering (like PBE), in contrast with measurements. In 
order to stabilize the correct G-type AFM arrangement 
a larger value of a is required. But these larger por- 
tions of exact exchange lead to a band gap significantly 
larger than experiment. The 'optimum' choice is proba- 
bly a=0.15 for which HSE delivers an AFM-G insulating 
solution with a band gap of about 0.7-0.8 eV (depend- 
ing on the structural details). For larger a the computed 
band gap is exceedingly large: 1.8 and 2.8 for a=0.25^- 
and a=0.35, respectively. 

The tendency of beyond-DFT methods to overestimate 
the bandgap of LaTiOa was already reported in litera- 
ture, based on SIC (1.7 eV-±22) and other HSE^ (1.7 eV 
using a=0.25, in agreement with our data) studies, and 
attributed to dynamical effects not included at this level 
of theory^. Furthermore, HSE tends to overestimate 
the magnetic moment of about 30 %, again in analogy 
with previous beyond-DFT studies. 

The MH like character of the band gap is evident by 
comparing the PBE and HSE DOS given in Fig|BJ the 
inclusion of non-local exchange split the t2 g band near 
Ep, thus opening a MH band gap between occupied and 
unoccupied ti g subbands. As expected the band gap in- 
creases with increasing a. The presence of an isolated 
peak on top of the valence band, well separated from 
the states beneath has been also detected by X-ray pho- 
toemission spectroscopy (XPS) experiments^^. The CT 
gap, defined as the energy separation between the O 2p 
states and the upper t2 g Hubbard band is also a depen- 
dent, and its value for the 'optimum' 0.15 choice, 4.7 eV, 
is in excellent agreement with experiment, 4.5 eV— . 

Finally, we underline that HSE is able to stabilize 
the correct orbitally-ordered state manifested by a chess- 
board G-type arrangement of differently ordered t^g 
cigar-lobes. We will come back to this point in the next 
section. 



3. d 2 : LaVOz 

LaVOs is another challenging material for conventional 
DFT: it is a i|g AFM-C Mott insulator, but DFT finds 
an AFM-C metal. The C-type antiferromagnetic spin or- 
dering is stabilized by the JT induced bond length alter- 
nations in the ab plane which cause the G-type orderings 
of d yz and d zx orbitals^. The experimentally observed 
MH and CT gaps are 1.1 and 4.0 eV, respectively^. 

Regardless the fraction of non-local exchange, HSE 
correctly finds a AFM-C MH insulating ground state, 
in which the gap is open between the lower and the 
upper MH t2 g band, similarly to LaTiO,-? (in PBE the 
t2 g band crosses the Fermi level, see Fig. [5J. The best 
agreement with experiment is achieved for a=0.10-0.15 
for which HSE delivers satisfactory values for both the 
MH (ps 0.8-1.4 eV for a =0.10 and a=0.15, respectively, 
as summarized in Table IXIII|) and CT gaps (w 4.4-4.9 
eV for a=0.10 and a=0.15, respectively). Similarly to 
all other theoretical DFT and beyond-DFT approaches, 
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Experimental Structure 



Optimized Structure 




-6-4-2 2 4 

Energy (eV) 



6 -4 -2 2 4 6 
Energy (eV) 



FIG. 8. (Color online) /-projected DOS of AFM-C ordered 
LaV03 with experimental (left) and relaxed (right) structure 
based on PBE and HSE (HSE-35, HSE-25, HSE-15, HSE-10) 
functionals. The shadow area indicates the total DOS. 
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FIG. 9. Bandstructure of LaVOa computed at HSE level 
(a=0.15) using the optimized structure. 



TABLE XIII. The band gap A (eV), magnetic moment m 
(/xs/V), magnetic energy (given with respect to the FM en- 
ergy, in meV) of LaVOa. calculated by PBE and HSE (HSE- 
35, HSE-25, HSE-15, HSE-10) using both the experimental 
and relaxed structures (Table |IV)) . Other theoretical values 
are also listed for comparison, along with the experimental 
measurements. 

Theory 
Optimized Structure 





HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 


A 


3.42 


2.43 


1.455 


0.885 


0.000 


m 


1.876 


1.855 


1.819 


1.782 


1.625 


A-AFM 


-73 


-54 


23 


43 


-77 


C-AFM 


-124 


-114 


-144 


-177 


-216 


G-AFM 


-96 


-98 


-30 


33 


137 






Experimental Struct 


ure 






HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 


A 


3.675 


2.535 


1.380 


0.810 


0.000 


m 


1.882 


1.858 


1.813 


1.774 


1.629 


A-AFM 


-2 


33 


16 


11 


-64 


C-AFM 


-105 


-119 


-151 


-179 


-124 


G-AFM 


-89 


-80 


-52 


-11 


203 






Other works 








LDA 


LDA+U 


GW 


HF 




A 


0.1" 


0.7 \ 0.92 c 
1.2 d 


2.48 c 


3.3 e , 0.9' 




m 


1.47 a , 1.85 6 


1.98\ 1.79 c 
1.70 d 


1.79 c 


1.8 e ,1.64 / 




A-AFM 


9° 


3.7 d 








C-AFM 


-35 a 


-38.3 d 








G-AFM 


17 a 


-U.8 d 









Experiment 


A 




1.1 s 


777 




1.3 h 


a RefM, 


6 Ref4£, c Ref.-ii, 


d Ref.^2, e Ref.^a, ^Refi^ 9 Ref.^, 


h Refi 






A- 


d 3 : LaCrOz 



HSE tends to overestimate the magnetic moment. It has 
been proposed that the origin of this discrepancy could 
by an unquenched orbital magnetization or spin-orbit in- 
duced magnetic canting^. 

The bandstructure of LaVC>3 computed for the repre- 
sentative case a = 0.15 is displayed in Fig|9l Also in 
this HSE is able to stabilize the correct G-type orbitally 
ordered state. This will be discussed in more details in 
the next section. 



Under equilibrium conditions P„ ma -distorted LaCr03 
exhibits a G-type AFM insulating ground state with the 
Cr +3 cation in the d 3 electron configuration t^g ■ The op- 
tical experiments by Arima et al. reported a coexistence 
of CT and MH like excitations in LaCr0 3 at 3.4 eVi 
These findings have been explained by several theoretical 
HF^S, LDA+US, GWli studies in terms of a significant 
mixing between Cr t2 g and O p states at the top of the 
valence band. In particular, the LDA+U study of Yang 
and coworkers has shown that the CT/MH character of 
the band gap is strongly U dependent: for small values 
of U (U<5 eV) the top of the valence band is mainly 
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Experimental Structure 



Optimized Structure 




-6-4-2 2 4 

Energy (eV) 



4 -2 2 4 
Energy (eV) 



FIG. 10. (Color online) Z-projected DOS of AFM-G ordered 
LaCrOs with experimental (left) and relaxed (right) structure 
based on PBE and HSE (HSE-35, HSE-25, HSE-15, HSE-10) 
functionals. The shadow area indicates the total DOS. 
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FIG. 11. Bandstructure of LaCrOs computed at HSE level 
(a=0.15) using the optimized structure. 



formed by t2 g states and the gap is predominantly MH, 
but for larger U (U>5 eV) the O p bands is progressively 
shifted towards higher energy thus reducing the size of 
the charge-transfer gap which become indistinguishable 
from the MH one. Our HSE calculations confirm this pic- 
ture as shown in the DOS plotted in FigfTOl the Op-Cr^ 
mixing at the top of the valence band increases with in- 
creasing a. As expected, a also influences the predicted 
band gap size which is found to be much smaller than 
experiment at purely PBE level (1.2 eV) and reaches the 
value 3.0 eV for cv=0.15, in good agreement with the 
reported optical gap. For larger a the gap starts to de- 



TABLE XIV. The band gap A (eV), magnetic moment m 
(fiB /Cr), magnetic energy (given with respect to the FM en- 
ergy, in meV) of LaCrOs. calculated by PBE and HSE (HSE- 
35, HSE-25, HSE-15, HSE-10) using both the experimental 
and relaxed structures (Table |V|. Other theoretical values 
are also listed for comparison, along with the experimental 
measurements. 

Theory 
Optimized Structure 





HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 


A 


5.475 


4.230 


3.000 


2.415 


1.245 


m 


2.866 


2.836 


2.790 


2.756 


2.643 


A A TH "TV T 

A-AFM 


-79 


-91 


-108 


-121 


-166 


C-AFM 


-160 


-184 


-221 


-245 


-309 


G-AFM 


-226 


-258 


-305 


-338 


-432 






Experimental Struct 


ire 






HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 


A 


5.460 


4.245 


3.030 


2.430 


1.245 


m 


2.868 


2.835 


2.784 


2.748 


2.626 


A-AFM 


-76 


-91 


-113 


-128 


-171 


C-AFM 


-170 


-203 


-249 


-281 


-375 


G-AFM 


-233 


-275 


-335 


-376 


-494 






Other works 








LDA 


LDA+U 


GW 


HF 




A 


1.40/3.4 e 


1.04 / , 1.40 e 


3.28* 


4.5" 




m 


2.56 e 


2.58 / , 3.00" 


2.38 J 


3.0" 




Experiment 


A 






3.4 C 






m 




2.45 a , 2.8 


b , 2.49 d , 2.63 





"RefiS 
"Ref.-SS. 



>Ref.- 



c Ref.-i, d Ref.J 



s Reft 



68 / 



Ref.-J 



l Ref.-i 



viate substantially from the measure value, and become 
exceedingly large for a=0.35 (see Table |XTV|) . The band- 
structure corresponding to the 'optimum' a=0.15 choice 
is displayed in Fig. [TTJ The G-type spin ordering is 
very robust at any value of a and the magnetic moment 
changes by only 0.2 \xb going from a—0 (« 2.6 (1b) to 
a=0.35 (~ 2.8 /is). Also in this case the electronic and 
magnetic properties obtained from the optimized struc- 
ture are essentially identical to those corresponding to 
the experimental structure. 

A different interpretation of the bandstructure and op- 
tical properties of LaCrOs was proposed in 2008 by Ong 
and coworkers who suggested that LaCrOs should not 
be considered a strongly correlated material^ These au- 
thors have attributed the 3.4 eV CT gap as the excitation 
from the top of the wide O p band below the t^g states 
to the bottom of the Cr d unoccupied band, and called 
for a new optical experiment to confirm the presence of 
a smaller MH gap of 2.38 eV open between Cr t^g and 
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Experimental Structure 



Optimized Structure 




Energy (eV) 



Energy (eV) 



FIG. 12. (Color online) /-projected DOS of AFM-G or- 
dered LaMnOa with experimental (left) and relaxed (right) 
structure based on PBE and HSE (HSE-35, HSE-25, HSE- 
15, HSE-10) functional. The shadow area indicates the total 
DOS. 
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TABLE XV. The band gap A (eV), magnetic moment m 
(/js/Md), magnetic energy (given with respect to the FM 
energy, in meV) of LaMnOs. calculated by PBE and HSE 
(HSE-35, HSE-25, HSE-15, HSE-10) using both the experi- 
mental and relaxed structures (Table IVI|I . Other theoretical 
values are also listed for comparison, along with the experi- 
mental measurements. 

Theory 
Optimized Structure 







HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 




A 


3.41 


2.47 


1.63 


0.75 


0.00 




rn 


3.78 


3.74 


3.67 


3.65 


3.52 


A- 


AFM 


-7 


-8 


-24 


3 


171 


C- 


AFM 


156 


182 


198 


368 


564 


G 


AFM 


161 


192 


208 


428 


899 








Experimental Structure 








HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 




A 


3.30 


2.40 


1.52 


1.10 


0.23 




in 


3.78 


3.73 


3.67 


3.62 


3.50 


A- 


AFM 


-4 


-11 


-28 


-44 


-63 


C- 


AFM 


164 


182 


198 


202 


209 


G 


AFM 


175 


195 


212 


216 


228 










Other works 








GGA 


GGA+U 


B3LYP 


HF 


GW 




A 


0.70" 


1.18" 


2.30 6 


3.0 r 


1.6 C , 1.68 d 




in 


3.33 a 


3.46 a 


3.80 6 


3.9 r 


3.16 c 






3.39 e 




3.77 7 


3.96' 


3.51 d 


Experiment 




A 




1.1 s , 


1.7 h , 1.9* 


, 2.0 ; ' m 






m 




3.87", 3.7°, 3.42 p 





a Ref.-^ 6 
Ref^21. 
Ref.-i^, 



'lief. 



c R e fIi, d RefM, e Ref.-^ 
Ref^, 1 Ref^, m Ref^ 
Ref^. 



., » RefiS, ° 



9 Ref.I, h 
Ref^I, p 



FIG. 13. Bandstructure of LaMnOs computed at HSE level 
(a=0.15) using the optimized structure. 



Cr e g bands, which would justify the green-light color of 
LaCrOs. We are not aware of more recent experimental 
data in support of this interpretation. 



5. d 4 .- LaMnOa, 

LaMnC>3 is one of the most studied perovskitc. Its 
properties have been widely studied both experimen- 
tally and theoretically as mentioned in the introduction. 
The initial tentative assignment of Arima and cowork- 
ers on the CT electronic nature of LaMnOs was succes- 



sively disproved and nowadays it is widely accepted that 
LaMnC>3 represents the prototypical example of a JT- 
distorted MH orbitally-ordered antiferromagnetic (type- 
A) insulato r 35 ' 160 ' 161 . In discussing the structural prop- 
erties we have underlined that LaMnC>3 is a very crit- 
ical case for conventional band theory due to the small 
but crucial JT distortions which are only marginally cap- 
tured by PBE. The drawbacks of standard DFT are also 
reflected in the electronic and magnetic properties sum- 
marized in Figfl2] and Table IXV1 especially for the theo- 
retically relaxed structure. Using the optimized geometry 
PBE favors the wrong magnetic ordering (FM) and stabi- 
lizes a metallic solution, whereas by adopting the exper- 
imental structure the correct AFM-A insulating ground 
state is stabilized, but the value of the band gap, 0.23 cV, 
is significantly smaller than the experimental one, 1.1-2.0 
eV (this is in agreement with previous studies^iSi) . This 
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Experimental Structure 



Optimized Structure 




-2 2 4 

Energy (eV) 



6 -4 -2 2 4 6 
Energy (eV) 



FIG. 14. (Color online) ^-projected DOS of AFM-G ordered 
LaFeC>3 with experimental (left) and relaxed (right) structure 
based on PBE and HSE (HSE-35, HSE-25, HSE-15, HSE-10) 
functionals. The shadow area indicates the total DOS. 



indicate that the JT distortions alone are sufficient to 
open up a band gap in LaMn03, but in order to predict 
a more accurate value it is necessary to go beyond DFT. 
In fact, turning to HSE the situation improves signifi- 
cantly and the results achieved within the theoretically 
optimized geometrical setup are essentially identical to 
those obtained for the experimental structure. The only 
significant difference regards the relative stability of the 
AFM-A ordering with respect to the FM one. For a=0.10 
the FM ordering is still more favored over the AFM-A one 
using the optimized geometry, but by adopting the exper- 
imental the AFM-A arrangement become the most stable 
one. For larger values of a both structural setups lead 
to essentially the same relative stability among all con- 
sidered spin arrangements. As expected, the band gap 
increases linearly with increasing mixing parameter and 
the best agreement with the measured values is reached 
again for a=0.15 (« 1.6 eV, well within the experimen- 
tal range of variation). The band gap is open between 
occupied and unoccupied Mn e g states which are almost 
completely separated from the other bands, as clarified 
in the bandstructure plot provided in Fig. 1131 The as- 
sociated orbitally ordered state will be presented in the 
next section. The HSE prediction for the Mn magnetic 
moment is in good agreement with low temperature mea- 
surements, 3.7-3.87 /tg i 137 ' 167 , and previous B3LYP data 
(~ 3.8 hb)2&^- We observe a small increase of the mag- 
netic moment with increasing mixing parameter, a gen- 
eral tendency noticed for the other LaMC>3 compounds. 
A more extensive discussion of the ground state proper- 
ties of LaMnOs can be found in our previous work o 35 ' 82 . 



TABLE XVI. The band gap A (eV), magnetic moment m 
(/Us/Fe), magnetic energy (given with respect to the FM en- 
ergy, in meV) of LaFe0 3 . calculated by PBE and HSE (HSE- 
35, HSE-25, HSE-15, HSE-10) using both the experimental 
and relaxed structures (Table IVTlj) . Other theoretical values 
are also listed for comparison, along with the experimental 
measurements. 

Theory 
Optimized Structure 





HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 


A 


4.680 


3.570 


2.460 


1.875 


0.660 


m 


4.198 


4.110 


4.001 


3.933 


3.719 


A- A 1 1V1 


-zoy 


-323 


-417 


A C7 
-4o f 


- 1 o 


C-AFM 


-530 


-653 


-832 


-947 


-278 


G-AFM 


-760 


-930 


-1166 


-1316 


-696 






Experimental Struct 


ure 






HSE-35 


HSE-25 


HSE-15 


HSE-10 


PBE 


A 


4.665 


3.570 


2.445 


1.875 


0.615 


m 


4.202 


4.111 


3.998 


3.927 


3.708 


A-AFM 


-251 


-321 


-427 


-511 


-9 


C-AFM 


-518 


-655 


-854 


-993 


-134 


G-AFM 


-742 


-930 


-1194 


-1372 


-552 






Other works 








LDA 


LDA+U 


GW 


HF 




Band gap 


0.0 a 


0.10\ 2.1° 


1.78 6 


4.0 9 




m 


3.5 a 


3.54 6 , 4.1" 


3.37 6 


4.6 9 





Experiment 


A 


2.1 c , 2.4 d 


m 


3.9 e , 4.6 / 



"Ref.S, 6 RefA "Ref.- 7 -, d Ref.i^, e Ref.-±^, / Ref.-^ 7 -, 9 Ref.^. 
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FIG. 15. Bandstructure of LaFeOs computed at HSE level 
(a=0.15) using the optimized structure. 
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d b : LaFeOj, 



The electronic configuration of Fe 3+ ion in LaFe03 is 



Experimental Structure 



Optimized Structure 



the high spin state (t 2g 



)(eV'). Below the rather high 



magnetic ordering temperature Tjv = 750 Kj^£& LaFcOs 
displays a G-type AFM spin ordering, and the d 5 spin 
saturation prevents the formation of orbital ordering. 
Arima 7 " reported that the spectrum of LaFeC>3 is sim- 
ilar to that of LaMn03, except for an increase of the 
insulating gap which is found to be 2.1-2.4 eV, about 0.5 
eV larger than the LaMn03 energy gap. The bandgap 
is opened between the predominantly O-p and Fe-e g va- 
lence band maxima and the lowest unoccupied tig band 
as shown in the density of states of Fig. [14] As such 
LaFe03 should be considered an intermediate CT/MH 
insulator, as originally suggested by Arima, who found 
almost identical CT and MH gaps 7 -. PBE does an ap- 
preciable job in predicting the correct AFM-G insulating 
ground state, though the value of the band gap, ~ 0.6 
eV, is significantly underestimated with respect to exper- 
iment (see the collection of electronic and magnetic data 
in Table IXVlj) . Similarly, the PBE estimates of the mag- 
netic moment, 3.7 \ib, is below the observed value. How- 
ever, it should be noted that the available low temper- 
ature experimental measures of the magnetic moments 
are very different, 3.9 ^et^ and 4.6 het^, thus a firm 
comparison is presently out of reach. 

The best agreement with the experimental gap is ob- 
tained also in this case for a=0.15 for which HSE gives 
a gap of about 2.4 eV, for both the optimized and ex- 
perimental structure (this is not surprising considering 
that in LaFcOs the optimized structure differences by 
less than 1% from the experimental one, as discussed 
previously). For this value of the mixing parameter we 
achieve an excellent comparison with photoemission data 
of Wadati et alr^, in terms of the position and charac- 
ter of the main peaks at -0.5 eV (Fe-e g , O-p), -2 eV (Fe- 
hg-O-p) and -6 eV (Fe-e g , O-p). These findings agree 
with the GW spectra computed by NoharaZi. By in- 
creasing the fraction of HF exchange the position of the 
lowest occupied ti g and e g states are gradually pushed 
down in energy and become progressively more localized 
whereas the position and bandwidth of the O-p band re- 
mains essentially unaffected. This leads to a worsening 
of the comparison with the experiment for a > 0.25. The 
a = 0.15 bandstructure is shown in Fig |15l Finally, we 
note that the energy separation between the unoccupied 
tig and e g states (the two lowest conduction bands, re- 
spectively, as indicated in Fig |T4]) . about 1.3 eV, is al- 
most independent from a and in good agreement with 
x-ray absorption spectroscopy^ and the GW— results. 



7. d b : LaCo0 3 

The complex magnetic behavior of LaCoOs represents 
a great challenge for theory. At low temperature LaCoOs 
is a diamagnctic insulator in which the Co 3+ arc aligned 
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FIG. 16. (Color online) /-projected DOS of non-magnetic 
LaCoOs with experimental (left) and relaxed (right) structure 
based on PBE and HSE (HSE-35, HSE-25, HSE-15, HSE- 10) 
functionals. The shadow area indicates the total DOS. 



TABLE XVII. The band gap A (eV) of non-magnetic 
LaCoOs . calculated by PBE and HSE (HSE-35, HSE-25, 
HSE-15, HSE-10) using both the experimental and relaxed 
structures (Table IVIII[) . Other theoretical values are also 
listed for comparison, along with the experimental measure- 
ments. 
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in the low spin (LS) state (^t )( e g)j with a total spin 
S=0. At about 100 K it undergoes a transition towards a 
paramagnetic state associated with magnetic excitations 
involving high-spin {{t2 g ){e 2 g ), S=2) and intermediate 
spin {{t\ g ){e g ), S=l) configurations, and at high temper- 
ature (T w 500 K) it shows a second magnetic anomaly 
associated with an insulator-to-metal transition ^ 172 ' 173 . 
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Experimental Structure 



Optimized Structure 
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FIG. 17. Bandstructure of LaCoOs computed at HSE level 
(a=0.05) using the optimized structure. 
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FIG. 18. (Color online) /-projected DOS of FM LaNi0 3 
with experimental (left) and relaxed (right) structure based 
on PBE and HSE (HSE-35, HSE-25, HSE-15, HSE-10) func- 
tional. The shadow area indicates the total DOS. 



These issues have been widely discuss in literature but a 
general consensus is still missing and their detailed un- 
derstanding remain highly controversia l 74 ' 77 . 

Standard LDA (or GGA) predicts a metallic and mag- 
netic ground opposite to experimen t 52 1 6 1 1 75 . Conversely, 
DFT+U can reproduce the correct non-magnetic insulat- 
ing state, but the results depends critically on the choice 
of U and the results seems to be strongly dependent on 
the specific computational schemes adopted. Small val- 
ues of U (< 2 cV) lead to the erroneous DFT-like solu- 
tion. It has been shown that the correct LS insulating 
solution can be obtained using rather different U, ranging 
from U«3 eV-S&S to XJ^8&Z&21. Hsu and coworkers^ 
have recently performed an optimization of the value of 
U based on an accurate account of the structural proper- 
ties, and show that the best agreement with experiment is 
achieved for a rather large U«8.2 eV— . A similar value of 
U has been also found by Laref et al. throughout the in- 
verse response matrices 7 - 7 -. Finally, using the unscreened 
hybrid functional PBEO scheme with the standard choice 
of the mixing parameter (0.25) Gryaznov et al. were able 
to find the correct LS state with a bandgap of 2.5 eV— . 
Our HSE results for a = 0.25 delivers a LS gap of 2.4 
cV, in excellent agreement with these PBEO results. 

In Table IXVlIl wc collect the values of the bandgap for 
the more stable S=0 HSE solution along with available 
experimental and other theoretical estimations. The best 
comparison with experiment is achieved for a rather small 
a = 0.05 for which HSE delivers a bandgap of about 
0.1 eV, in good agreement with optical measurements, 
0.1-0.3 eV— i 1 - 7 - 4 - (photoemission data of Chainani et al*2^ 
give a somehow larger gap of about 0.6 eV). We remind 
that this value of a leads to the most accurate optimized 
geometry, as discussed in the previous section (see Fig. 
|2"2"|) . From the density of states shown in FigJTJ5]we evince 



that the gap is opened between valence band mixed O p 
and Co d states and empty d-like Co states, in agreement 
with the DFT+U and PBEO results mentioned above. 
The effect of the inclusion of a fraction of HF exchange 
is the splitting of the occupied ti g manifold and the e g 
states (this forms a continuous band which crosses the 
Fermi energy at PBE level). The valence band DOS is 
characterized by three main regions located at -1 eV, -3 
cV and -5 cV, reproducing well the XPS±2£ and GWli 
spectra. Equally satisfactory is the distribution of the 
conduction band states, with the Co and La d states 
centered at sa 2 eV and 7-9 eV, respectively. 

The bandstructure plotted in Fig[T7] provides further 
evidence for the large degree of hybridization of the top 
of the valence band and the rather dispersive character 
of the lowest e g unoccupied states. On the basis of this 
analysis LaCoOa can thus be considered to be predomi- 
nantly a CT like (O p — > Co d) insulator (in agreement 
with the initial assignment by Arima 7 -) but a Mott mech- 
anism is necessary to split apart the Co d band crossing 
Ep at PBE level, possibly indicating minor t^g — > e g MH 
like excitations, which have not been specifically investi- 
gated so far, by both theory and experiment. 



8. d 7 : LaNiOi 

LaNiOa is a weakly correlated PM metal in which 
the Ni +3 ion possesses the low-spin 3d configuration 
(tig The electron-electron correlation associ- 

ated with the partially filled Ni 3d 7 shell is inhibited by an 
efficient electrostatic screening, originated by the strong 
Ni 3d-0 2p hybridization (relatively small Ni-0 distance), 
and d-d hybridization (large valence rf-bandwidth) 7 - ' 47 ' 79 . 



21 



TABLE XVIII. The band gap A (eV) and magnetic mo- 
ment m (hb /Ni) of FM ordered LaNiC>3 , calculated by PBE 
and HSE (HSE-35, HSE-25, HSE-15, HSE-10) using both the 
experimental and relaxed structures (Table [Kj. Other theo- 
retical values are also listed for comparison, along with the 
experimental measurements. 
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The electronic structure of LaNi03 has been recently 
extensively investigated and thoughtfully discussed by 
the group of J. M. Rondinelli 7 -^ using an array of sev- 
eral above-standard first principles methods including 
LDA+U, PBE0 and HSE. 

Our HSE results (summarized in FigQ~8] and Table 
IXVIIip reproduce the trends observed by Rondinclli and 
coworkers and support their conclusions: 

(i) Conventional DFT works fairly well as it provides 
a correct non-magnetic metallic solution. This was al- 
ready pointed out in precedent work o 176 ' 177 . We should 
however note that for this specific case structural effects 
are extremely important in the determination of the rel- 
ative stability between the non-magnetic ground state 
and the competing FM solution. Using the experimental 
structure PBE favors the non-magnetic case by about 
130 meV/f.u., but adopting the PBE-optimized struc- 
ture the FM ordering become the most stable solution 
by about 110 meV/f.u.. This should be attributable to 
the PBE overestimation of the volume (+2.3 %), as all 
other structural properties are described by PBE with 
an error smaller than 1% (see Table HX|) . The compari- 
son with PES data gives further support to the quality 
of the DFT performance, as discussed in Refl79l. Minor 
differences have been observed between LSDA and PBE, 
relative to the width of the valence band which is better 
described at LSDA level. The bandstructure computed 
at PBE level for the most stable non-magnetic solution 




FIG. 19. Bandstructure of LaNiOs computed at PBE level 
using the optimized structure. 



given in Fig|l~9lshow evident similarities with the LaCo03 
bands. The major difference is the downward shift of the 
e g manifold at the bottom of the conduction band which 
now crosses the Fermi level and get mixed with the lower 
laying occupied t2 g orbitals. The strong d-p and d-d hy- 
bridization is reflected in the highly dispersive character 
of the valence bands, in accordance with the DOS. 

(ii) HSE, similarly to PBE0 and DFT+U 7 ^, delivers a 
very lacking picture: LaNiOa is described as a FM half- 
metal with a magnetic moment m of about 1 /is. m 
increases gradually as a function of a and reaches the 
value 1.3 /is for cv=0.35. The deficient HSE results is 
mostly due to an excessive downward shift of the ti g man- 
ifold (this increases the overall bandwidth with respect 
to PBE), a strong depletion of Ni d states on top of the 
valence band, and a much too large exchange splitting. 
Clearly, the relative strength of these effects increases 
with increasing a, as clarified in FigdH] This makes the 
comparison with the PES data much worse, in terms of 
both bandwidth, and number and positions of the main 
peaksi 7 ^. The fundamental failure of hybrid functional 
for itinerant magnets was already reported by Paier et 
al. for bulk Fc, Co and Nii22. 

Though conventional DFT leads to a decent account 
of the ground state of LaNiOa it should be emphasized 
that LaNiOs is experimentally recognized as being a cor- 
related metals, with important dynamical correlation ef- 
fects associated with the Ni e g orbital*^ which cannot 
be captured at DFT level. More suitable methodolo- 
gies such as dynamical mean-filed theory are need to ap- 
preciate the fundamental nature of LaNiOs, as recently 
dcmonstratecU 7 ^^. 
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Experimental Structure 



Optimized Structure 
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FIG. 20. (Color online) /-projected DOS of non-magnetic 
LaCuOs with experimental (left) and relaxed (right) structure 
based on PBE and HSE (HSE-25, HSE-15, HSE-10) function- 
als. The shadow area indicates the total DOS. 
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FIG. 21. Bandstructure of LaCuOs computed at PBE level 
using the optimized structure. 



9. d s : LaCuOa 

LaCuC>3 is a PM metali. Cu 3+ ions are formally in 
the low spin configuration (t^g )(e^) (the t2 g shell is 
fully occupied and the e g orbitals are half- filled) , but it 
has been argued that this d 8 state is strongly hybridized 
with the d 9 L configuration in which L denotes a lig- 
and hole, thus suggesting the existence of orbital de- 
generacy associated with significant Cu-Cu many-body 
excitation3i2i~-i£2. This is another challenging case both 
for theory (orbital degeneracy and dynamical correla- 
tion) and experiment (it is very difficult to to prepare a 



TABLE XIX. The band gap A (eV) of non-magnetic 
LaCuOs, calculated by PBE and HSE (HSE-35, HSE-25, 
HSE-15, HSE-10) using both the experimental and relaxed 
structures (Table |X)) . HSE-35 favors an FM-ordered ground 
state with m—1.197fiB , similarly to the LDA+U calculation 
of Ref. [4r3 . Other theoretical values are also listed for com- 
parison, along with the experimental measurements. 
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stoichiometric tetragonal phase of LaCuOa without oxy- 
gen vacancies). Thus, the final methodological comments 
given for LaNiC>3 on the necessity to employ many-body 
schemes to achieve a fundamentally more accurate theo- 
retical description remain valid for LaCuOs as well. 

Our PBE and HSE results are collected in Table IXIXl 
and Figs. [50] and [5TJ In agreement with the results of 
M.T. Czyzyk and G.A. Sawatzky^ we find that stan- 
dard DFT finds the correct metallic non-magnetic ground 
state. The DOS (Fig. [2H|) is characterized by a wide 
band crossing the Fermi level formed by Cu d (primar- 
ily e g ) and O p states, associated with an highly dis- 
persive bands (Fig. [2T]) . In analogy with H F 50 ' 181 and 
LDA+U^ calculations also within HSE the G-type AFM 
insulating state is lower in energy than the non magnetic 
metallic state, in contradiction with experiment. Here 
we only report the results for the non-magnetic solution. 
From the DOS shown in FigJ^Ol we infer that the elec- 
tronic structure stays almost unchanged with respect to 
the non-magnetic PBE case. The only notable difference 
is a progressive downward shift of the t2 g Cu states with 
increasing a and a gradual further broadening of the Cu 
d/O p band crossing the Ep. 

We conclude this section by providing in Fig. [22] 
a schematic graphical interpctation of the comparison 
between computed and measured structural, electronic 
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and magnetic properties, given in terms of the obtained 
MARE. A more elaborated discussion will be developed 
in the next section. 

IV. DISCUSSION 

From the analysis of the structural, electronic and 
magnetic properties developed in the previous section we 
have derived a set of 'optimum' values for the mixing pa- 
rameter (a^f t E ) for which HSE (and in two cases PBE, i.e. 
a=0) delivers a substantially correct and quantitatively 
satisfying description of the LaMC>3 series (within the 
limits discussed previously). This set of oS E parame- 
ters, collected in Table IXXl includes a=0.25 (for the wide 
band gap insulator LaScOa), a=0.1-0.15 (for the MH 
and intermediate MH/CT insulators LaTiOs-LaFcOa), 
a=0.05 (for the small band gap CT insulator LaCoOs), 
and a=Q (for metallic LaNiOa and LaCuOs). As already 
reported in Sec lIII B"2| it is important to underline that 
for LaTiC>3 the overall best quantitative agreement with 
experiment is achieved for a=0.l (the computed band 
gap is ~ 0.2 eV, almost identical to the measured value), 
but the incorporation of this fraction of exact exchange 
in HSE leads to the wrong magnetic ordering, AFM-A 
instead of AFM-G. 

It is instructive at this point to compare the set of pa- 
rameters aop S t E with the optimum values of a obtained 
throughout the inverse dielectric relation a e op ° t ~ j— in- 
troduced in the computational section (Eq. [7]) and de- 
rived by mapping hybrid DFT with GW. The measured 
dielectric constant £ocr^ and the corresponding — val- 
ues are also listed in Table IXXl The nice correlation be- 
tween these two sets of numbers can be appreciated in 
Fig. [23] where we display dJ) p s t E and c\£x. The two curves 
follow a very similar behavior characterized by an initial 
large value of a for the poorly screened d° band insula- 
tor LaSc03 followed by a plateau of similar values in the 
range d 1 (LaTiOs) — > d 5 (LaFeOa) and finally a sharp de- 
crease towards the more strongly screened metallic com- 
pounds characterized by a completely filled t2 9 manifold. 
For LaNiC>3 and LaCu03 the optimum value of a is zero 
(not shown). Thus, the a pt curve derived from the HSE 
fitting procedure conducted by computing a wide set of 
structural (volume, cell shape, JT distortions, atomic po- 
sitions) electronic (band gap and DOS) and magnetic 
(spin ordering, magnetic moment) properties as a func- 
tion of a reflects well the evolution of the screening prop- 
erties across the LaMC>3 series represented by the dielec- 
tric function f x . However, from a quantitative point of 
view the two sets of value differ by about 0.07, as clarified 
graphically by the open squares in Fig. [531 in order to 
achieve a good quantitative match between the cc^| E and 
a^ t curves it is necessary to shift downwards the latter 
by about 0.07. This behavior is attributable to two main 
reasons: (i) The relation a e op ° t « ^-) identifies a propor- 
tionality between a e ^ t and j— , not an exact equality (the 
factor of proportionality is not exactly 1, as discussed in 



Refs. Ill7l4ll9l) ; (ii) As already mentioned before, Eq. [7] 
holds for standard unscreened hybrid functionals such as 
PBE0. HSE is a range-separated screened hybrid func- 
tional which contains already a certain degree of screen- 
ing (controlled by the screening factor /i). Therefore it is 
expected that the optimum a derived for PBE0 (a^ E0 ) 
will be smaller than the corresponding /^-dependent HSE 
value (a^p t E )^i^. Needless to say, that in absence of a 
systematic study of the role of /i it is very difficult to 
quantify its effect on a^ s , E . We leave this issue open for 
future studies. 

Recently, Alkauskas et al. has proposed that an op- 
timal mixing coefficient can generally be found for any 
material using the formula*!^: 

^Expt /\ semilocal 

<t = 1 (8) 

where A Expt and A scmllocal indicate the experimen- 
tal and semilocal (GGA/LDA) band gap, and k = 
dA(a)/da (A(a) represents the variation of the band 
gap as a function of a)ii£. This practical relation takes 
advantage of the linear relation between A and a, which 
holds true as long as the valence band maxima and con- 
duction band minimum (and their associated wave func- 
tions) do not change much by changing a^^. In practice, 
if the experimental bang gap is known, it is sufficient to 
perform only one hybrid functional calculation for an ar- 
bitrary value of a and one semilocal calculation to derive 
the optimum value of a. Clearly, this empirical proce- 
dure does not guarantee that the so obtained optimum 
a is the best choice for what concerns the other ground 
state properties. We have already addressed this issue 
for LaTiOs for which a=0.1 give the best band gap but 
leads to the incorrect magnetic ordering. 

The changes of the bang gap as a function of a for the 
series LaScOvLaCoOs are reported in Fig. [M] The lin- 
earity is well preserved for most of the materials with the 
exception of the small band gap compounds LaTiO"3 and 
LaCoOs for which a sudden change of k = dA(a)/da 
is observed for a critical value of a. A departure from 
linearity is also found for the JT/MH insulator LaMn03 
if we adopt the fully relaxed structure (full line). This 
is due to the fact that the cooperative JT distortions, 
which contribute to the opening of the band gap, do not 
change linearly with a (as such, the associated wave func- 
tion will change with a). Indeed, by keeping the atomic 
coordinates fixed to the experimental positions the gap 
grows linearly by increasing a (dashed line). 

The values of a^ pt obtained from this set of curves are 
indicated with empty (red) circles and included in Table 
IXXl For some materials we provide more than one opti- 
mum mixing parameter since different experimental gaps 
are reported in literature, as reported in the previous sec- 
tion. Not surprisingly, we find that the values of a^ pt are 
very similar to the corresponding a^ s t E data, with the 
exception of LaTiOs and to a lesser extent LaMnOs, and 
correlates well with a t ( ^ tl as visualized in Fig. 1231 
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FIG. 22. (Color online) Summary of the MARE for the structural properties (top panel), band gap A (middle panel) and 
magnetic moment m (lower panel), at PBE and HSE level. For the band gap A and the magnetic moment m the MARE 
are indicated by the numbers associated to each bar. A few specifications for the labels 'OK' and 'wrong': (i) LaSc03, m: 
all methods correctly predicts a non-magnetic ground state; (ii) LaCo03, m: all methods correctly predicts a zero magnetic 
moment; (hi) LaNiC>3, A: PBE is the only approach which correctly find a metallic solution; (iv) LaNiC>3, m: all methods 
wrongly predict a magnetic ground state; (v) LaCuC>3, A: all methods correctly predicts a metallic solution; (vi) LaCuC>3, 
m: PBE and HSE (0.05, 0.15, and 0.25) correctly predicts a zero magnetic moment, whereas HSE 0.35 wrongly stabilized a 
magnetic ground state. 



TABLE XX. Comparison between the set of optimum mixing factors a for the 3d perovskite series LaMOs (M=Sc-Cu) 
computed throughout the HSE fitting procedure developed in Sec. IHII and those obtained using the relations a^JJJ = I/too (Eq. 

and OL^pt — — — ~~~~t (Eq. [8|). The values of the dielectric constant, given in brackets, are taken from Ref ll84l . 
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Now, with the rough-and-ready set of optimum HSE- 
fitted values a^ p s t E we conclude this paper by providing 
a general picture of the variation of the properties of 
LaM O3 from M=Sc to M=Cu, by comparing our com- 
puted results with the available experimental data. Fig- 
ureUHlshows the remarkably good agreement between the 
calculated and measured values of the Volume (V) , tilting 
angle (6), JT distortion, band gap (A) and magnetic mo- 
ment (m) . The correlation between V and Rm , as well as 
between 6 and t was already discussed at the beginning 
of Sec. IIII Al The variation of the magnetic moment as 
a function of M can be easily understood in terms of the 
progressive ti g and e g band filling in the high-spin com- 
pounds LaTi0 3 (t2g^, to=0.51 /jb), LaV0 3 (t^g^, m=1.3 
/i B ), LaCr0 3 (t 2s m , m=2.63 /i B ), LaMn0 3 (t 2g ttt e g t , 
to=3.66 /j, B ), and LaFe0 3 (t 2g m e g n , m=3.9-4.6 /jB). 
As already specified the experimental and computed 



magnetic moments of LaCo0 3 should be taken with a 
certain caution. LaNi0 3 and LaCu0 3 shown a non- 
magnetic behavior at PBE level only (we have already 
reported that the a small magnetic moment of 0.169 /jb 
for LaNi0 3 is found for the fully relaxed structure). 

The variation of the band gaps with the M species 
shown in Fig. is rather complex and in line with the 
earlier observation of Arimaii^: we observe a general 
tendency of the MH gap to increase as the TM atomic 
number increases, whereas the CT gaps follow an oppo- 
site behavior. As expected, there is an apparent corre- 
lation between the trend of the band gaps and the op- 
timum fraction of exact exchange displayed in Fig. [2^1 
especially a^ pt . In LDA+U language this behavior is in- 
terpreted as a correlation between the strength of the 
effective Coulomb interaction U and the band gap (this 
become particularly evident by comparing the A curve 
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FIG. 23. Graphical interpretation of the optimum values of a 
listed in Table IXXl showing the correlation between the HSE 
fitted parameters (HSE fit), the inverse dielectric constant 
relation (1/eoo), and a^ pt (Eq. [HJ. The light-gray squares 
represent the l/e^o values shifted by 0.07. This shift roughly 
reflects the amount of screening incorporated in HSE via the 
screening factor \i as compared to fully unscreened PBEO (see 
text). 

in FigJ55] with Fig. 2 in Ref. HH showing the changes of 
the effective U with respect to M). 

In Fig. (21]wc collect the bandstructures of LaMC>3 ob- 
tained using the optimum a^f t E values showing the vari- 
ation of the electronic dispersion across the whole series. 
Starting from the d° band insulator LaSc03 the addition 
of one d electron creates an highly localized t 2g state 
right below Ep in LaTiC>3. The gradual filling of this 
tig manifold leads to a continuous increase of the band 
width from t 2g 1 (LaTiOs) to t 2g 3 (LaCrOs), connected 
with a gradual increase of the crystal field splitting. In 
LaMnC>3 the fully occupied t 2g band is pushed down in 
energy and the valence band maxima is dominated by 
the half-filled e g l subbands. The e g orbital gets com- 
pletely filled in LaFeC>3 which is the last member of the 
series having a predominantly MH gap. The inclusion 
of one additional electron yields a sudden change of the 
bandstructurc manifested by a high increase of p-d hy- 
bridization and bandwidth around Ep. which finally lead 
to the onset of a metallic state in LaNiC>3 and LaCuC>3. 

Three members of the LaMC>3 family (LaTiC>3, 
LaVC>3, and LaMnOs) are known to display orbital- 
ordering (00) associated with the partially filled t 2g and 
e g orbitals located on top of the valence band (these 
states are demarcated by thick lines in Fig. 121)]) . A visual 
representation of the OO states derived from the 'opti- 
mum' HSE wavefunctions is shown in Fig |27l in terms 
of charge density isosurfaces of the highest occupied d 
states. In the following we describe briefly the most im- 
portant characteristics of the observed 00 states. 

(i) In LaTi03, where the 00 originates from the 
single t2 g electron, the lobes have a quasi cigar-like 
shape with asymmetric contributions along the two 
main directions, indicating an almost identical occu- 
pation of the three xy 7 xz, and yz t 2g shells. Co- 




mixing factor a 



FIG. 24. Change of the band gap as a function of a (think 
lines) and optimum values of a (circles) obtained throughout 

the practical formula a£~ pt =— — ——4 (Eq. ©■ 

planar lobes are arranged in a chessboard-like way with 
a sign alternation along z, in good agreement with pre- 
viously reported theoretica l 103 ' 152 ' 185 and experimental 
work o 131 ' 186 . There is a clear connection between this 
chessboard-like Ti d 1 ordering and the JT structural in- 
stability, which is manifested by the tendency of the oc- 
cupied t2 g state to lie along the longer Ti-0 bond. This 
also explains why the chessboard-like 00 in LaTi03 is 
not as much evident as in LaMnOs: in LaTi03 the dif- 
ference between the distinct Ti-0 bondlcngths Ti-O s Ti- 
O m Ti-O/, quantified by the JT parameters Q2 and Q3, 
is about one order of magnitude smaller than in LaMn03 
(see Table |m] and Table EJ). 

(ii) The V 3+ ions in LaV03 can accommodate two elec- 
trons in the three possible orbital states d xy , d xz , and 
d yz . The spins are arranged according to the C-type 
ordering, whereas the 00 state is found to be G-type, 
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FIG. 25. Trend of selected structural (Volume V, tilting an- 
gle 6, and JT distortions Q2 and Q3), electronic (bandgap 
A), and magnetic (magnetic moment m) quantities along the 
Lai?03 series from M=Sc to M=Cu. We also show the trend 
of the tolerance factor ^(Ra+RoVv^Rm+Ro), where Ra, 
Rm and Ro indicate the ionic radius for La, M=Sc-Cu and 
O, respectively, as well as Rm- For LaTi03 we used a=0.1. 
The character of the insulating gap is also indicated (BI = 
band insulator, CT = charge transfer, MH = Mott-Hubbard, 
CT/MH = mixed CT and MH character). 



in accordance with the Goodcnough-Kanamori rules^ 
and in agreement with x-ray diffraction^ 8 - and previous 
GGA^ and HE 9 ^ calculations. The distribution of the 
t2 g orbitals in the G-type 00 state follows the cooper- 
ative JT-induced V-0 bond-alternation in the xy plane 
and along the z axis, i.e. the tig charge density in one spe- 
cific V site is rotated by 90° with respect to that in the 6 
neighboring V sites (four in-plane and two in the adjacent 
vertical planes). As already observed for LaTiOa, the tig 
orbitals are preferentially occupied along the long-bond 
direction. 

(hi) The C-type OO in LaMnOa, originating from 
the singly occupied e g state of the Mn +3 3d elec- 
trons in the high-spin configuration t2 g 3 e g 1 has been 
extensively studied both experimentall y 164 ' 189 ' 190 , and 
theoreticall y 160 ' 161 . We have also recently addressed this 
issue throughout a maximally localized Wannicr func- 
tions representation of the e g states^. This C-type OO 
state can be written in the form \9) = cos||3z 2 — r 2 ) + 
sin||a; 2 — y 2 ) with the sign of 9 ~ 108° alternating along 
x and y and repeating along z, as correctly represented 
by our HSE charge density plots. 



For comparison we show in Fig. [28] the similar charge 
density isosurfaces calculated for LaCrOs and LaFeOs, 
in which the half-filling of the tig and e g orbitals inhibit 
any orbital flexibility. As expected there is no sign of 
orbital ordering. 

We conclude this section with a comparison of the cal- 
culated density of valence and conduction states with 
available x-ray photoemission and x-ray adsorption spec- 
tra. This is summarized in Fig. 1291 The overall pic- 
ture is satisfactory in terms of bandwidth, intensity and 
peaks position, though some sizable deviations are visi- 
ble for LaCrOs, LaFeOs and for the metallic compounds 
LaCuOs and LaNi03. These differences between theory 
and experiment should be attributable to the approxima- 
tions included in the adopted computational scheme, as 
mentioned in the previous section. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have studied the ground state proper- 
ties of the perovskite series LaM03 by means of screened 
hybrid DFT following the HSE formulation, based on 
the inclusion of a fraction of exact HF exchange in the 
short range Coulomb kernel of the underlying semilocal 
PBE exchange-correlation functional. In particular, we 
have investigated the role of the HSE mixing parameter 
a on the (i) structural parameters (Volume, JT/GFO 
distortions, lattice parameters and unit cell symmetry) 
(ii) electronic character (band gap, DOS, and bandstruc- 
ture), and (iii) spin orderings and magnetic moment, as 
a function of the gradual filling of the d manifold from 
LaSc0 3 to LaCu0 3 {d 8 : t 2g Q e g 2 ). 

On the basis of a computational fitting of the most 
relevant experimentally available data we have derived a 
set of mixing factors which leads to an accurate qualita- 
tive and quantitative description of the physical behav- 
ior of all members of the representative LaA/03 fam- 
ily. It is found that, apart from LaScOs, the 'opti- 



mum' values of a (a^ s , E ) are significantly smaller than 
the routinely used standard choice a=0.25, especially 
for the end members ( LaScOs: cv^f t E 



a"^r=0.25; LaTi0 3 



LaV0 3 : a 



,„HSE_ 
3- « op t - 

HSE_ 
opt 



0.10-0.15; 



LaFe0 3 : o£if=0.15; LaCo0 3 : o££ K =0.05; LaNi0 3 & 



LaCr0 3 , 

y HSE„ f 
opt 



LaMnOs, and 



LaCuOs: ay^. =0.0, i.e. for these two cases PBE is bet- 
ter than HSE ) . This can be understood by correlating the 
changes of a^f' from Sc to Cu with the corresponding 
values of the inverse dielectric constant 1 /coo , and by con- 
sidering that a certain degree of screening is already in- 
cluded by construction in the HSE functional throughout 
the screening length /1, at variance with the unscreened 
parent hybrid functional PBE0 (for which /i = 0). This 
suggests that the 'optimum' value of a in HSE should be 
smaller than the corresponding PBE0 one: in our spe- 
cific case it is proposed that the difference between a^ s ( E 
and aoj^ E0 should about 0.05-0.07, but a more detailed 
analysis on the influence of fi is required to achieve more 
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FIG. 26. (Color online) Summary of the HSE electronic dispersion relations showing the complete trend from the band insulator 
LaScOs to metallic LaCu03. The thik (red) lines demarcate the d bands responsible for the observed orbital-ordering in LaTiOs 
(tag), LaV03 (ta B ) and LaMnC-3 (e 9 ). 



accurate and comprehensive conclusions. 

An alternative way to obtain a set of optimum a is 
the fitting of the bandgap only, by applying the practical 
recipe represented by Eq. [5] However, though this proce- 
dure has the clear advantage of reducing considerably the 
computational cost, it can lead to an erroneous descrip- 
tion of other properties (for example the best-bandgap 
a in LaTiC>3 stabilizes the wrong spin ordering) and can 
only be applied under the assumption that the wavefunc- 
tion does not change with a. 

For what concerns the description of the modulation of 
the electronic and magnetic properties across the LaA/03 
series, we found that for all insulating compounds (M=Sc 
to Co) HSE is capable to capture the localization of the 
t eg /e g orbitals and to provide a consistent and quantita- 
tively satisfactory description of all considered quantities, 
thereby improving the deficient DFT based predictions. 

For the structural properties, on the other hand, 
PBE performs rather well, delivering optimized geome- 
try within 1%. The only exception are the JT parameters 



in LaMn03, which PBE finds 60% smaller than experi- 
ment. HSE cures this limitation, reproducing quite well 
the critical JT and GFO structural instabilities, and, in 
general, its application improves even further the overall 
agreement with experiment as compared to PBE. 

The complex nature of the PM correlated metals 
LaNiOa and LaCuOs, end members of the LaAfC>3 se- 
ries, is only marginally accounted for by PBE and rather 
poorly treated at HSE level. This is mostly due to un- 
derlying dynamical correlation effects which cannot be 
easily treated at DFT/HF level. For these compounds, 
PBE might be considered to be a good starting point for 
more elaborated many-body approaches. 
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FIG. 29. Comparison between experimental (blue squares) 
and calculated (red full lines) valence and conduction band 
spectra at the 'optimum' value of the a parameter: (i) 
LaSc0 3 : a = 0.25; (ii) LaTi0 3 : a = 0.15; (iii) LaV0 3 : 
a = 0.15; (iv) LaCr0 3 : a = 0.15; (v) LaMn0 3 : a = 0.15; 
(vi) LaFe0 3 : a = 0.15; (vii) LaCo0 3 : a = 0.05; (viii) 
LaNi0 3 : a = 0; (ix) LaCu0 3 : a = 0. The calculated and 
measured spectra have been aligned by overlapping the va- 
lence band maxima and conduction band minima. The ex- 
perimental data are taken from the collection of spectra pre- 
sented in Refl7ll. originally publis hed i n separated article s: (i) 
LaScOs: RefQ; (ii) LaTiQ 3 : Reffl9ll: (iii) LaVO s : Reffl92l: 
(iv) LaCrOs: Refill: (v) LaMn0 3 : Refill: (vi) LaFe0 3 : 
Reffml: (vii) LaCoOs: Refll28l: (viii) LaNi0 3 : Reffl9i: (ix) 
LaCuOs: Reffl8ll. 



